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SYNOPSIS 

"defect energetics in the alloy semiconductors" 

A Thesis Submitted 

in Partial Fulfilment of the Requirements 
For the degree of 

DOCTOR OF PHILOSOPHY 
by 

AMITA DAS 
to the 

Department of Physics 
Indian Institute of Technology - Kanpur 

January 12, 1990 

The alloy semiconductors are extensively employed in various 
electronic devices. They offer a possibility of engineering the 
material parameters of a semiconductor, notably the band gap and 
the electron mobility. 

Point defects however are invariably present in such 
semiconductors. These defects could be native; vacancy and 
anti-site; or foreign: donor, acceptor, transition metal atom etc. 
In some cases they are responsible for device degradation whereas 
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in others they are deliberately introduced to serve a specific 
obj ective . 

The main thrust of this thesis is the study of defect 
energetics in III-V alloy semiconductors. The alloy semiconductors 
(e.g. Ga A1 As or Si Ge ) are substitutionally disordered and 

1 X X X jL*~**X 

may be viewed as a generalization of the ordered III-V (e.g. GaAs) 
or group IV (e.g. Si) semiconductors. We concern ourselves with 
developing a formalism for defect energetics in such disordered 
hosts as well as conducting calculations for some typically 
important defect centers. It is possible to specialize the 
framework we develop to study defect energetics in ordered 
semiconductors. In a few instances we take advantage of this 
flexibility and calculate and comment on defect centers in Si and 
GaAs. Equally significant is the issue whether interpretation and 
analysis of experiments for ordered systems can be extended to 
alloy and other disordered semiconductors in an unmodified 
fashion. We address this issue in the context of Deep Level 
Transient Spectroscopy (DLTS). We also discuss a complementary 
technique : Paired Temperature Spectroscopy (PATS). 

As noted earliar, the host system Ga A1 As for instance, is 

1-X X 

substitutionally disordered with the cation site being randomly 
occupied by Ga or A1 atom. We employ the sophisticated Coherent 
Potential approximation (CPA) to describe the host electronic 
structure: We present examples of the CPA in a classical context. 
A new derivation of the Weiss molecular field result is presented^ 
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A deep level in the forbidden gap of a semiconductor may 
arise in the presence of a point defect. We present model single 
band Hamiltonians. We analyse the scalar version of the Koster 
Slater equation which serves to locate the deep level. Further, 
with the help of the Hellmann - Feynman theorem we calculate total 
energy changes due to the presence of point defects^. 

We develop an approximate version of CPA for the realistic 
multiband systems. This enables us to make quantitative 
predictions for defects in several alloy semiconductors and in 
particular we discuss a donor related level in Ga A1 As ( DX 

1-X X 
3 

center ) . In our calculations we have viewed this defect as a 
simple donor, which gives rise to a deep level on account of the 
disorder in the host structure. Our results compare well with the 
experiments . 

We have also studied helium like impurities (impurities 
having two excess valence electrons than the host) in the context 
of Effective Mass Theory (EMT). The chemical shift of the level 
has been obtained by introducing a cavity. Inside the cavity the 
excess electrons behave as in free atoms; outside, the 
semiconductor matrix dominates. We apply this formalism to a 
number of defects (i) chalcogen impurities in silicon (ii) DX 
center viewed as a negatively charged donor (iii) EL2 center in 
GaAs viewed as an anti-site defect. Our results are in 

4 - 

considerable agreement with the experiments . 

Finally, we have concentrated on the analysis of experiments 
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such as Deep Level Transient Spectroscopy (DLTS) and Paired 
Temperature Spectroscopy (PATS). These are used extensively for 
the determination of activation energies of the defect levels. We 
point out that the expression used to get the value of activation 
energy is valid only for ordered systems. This might result in a 
systematic discrepancy between the actual and the extracted value 

5 <5 

of the activation energy ' . We have also focussed on the effect 
of disorder on the M-center, a metastable defect in InP and its 
alloys^ . 
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CHAPTER 1 : INTRODUCTION 

1.1: Prelininary renarks 

Semiconductors play an important role in technology. There 
exists a great variety of them ranging from elemental to alloy and 
even amorphous semiconductors. Members of the group IV in the 
periodic table, like Ge and Si are semiconductors. The elemental 
semiconductors have been known for some time and are well studied. 

The compound semiconductors offer a wider choice of materials 
with differing characteristics such as band gap, lattice 
parameter, electron mobility etc. The popular compound 
semiconductors are the III - V or II - VI types. These are the 
compounds of elements of group III and V (for example GaAs ) and 
group II and VI (for example CdTe ) respectively. The relationship 
between the band gap and the lattice parameter of some of these 
compounds is presented in figure 1.1 . 

The alloys of the above mentioned elemental and compound 
semiconductors also have semiconducting properties and are known 
as alloy semiconductors. Similar to the compound semiconductors, 
the alloy semiconductors can also be classified on the basis of 
the periodic table: the III-V alloys (Ga^_^Alj^As , 
Ga ^In As, P ), the II-VI alloys (Hg,_ Cd Te) and the group IV 
alloys (Si Ge )- Another mode of characterising them is by the 

i-X X 



1-2 


presence of the number of constituents: binary (Si^j^Ge^^), ternary 
(Ga A1 As), and quaternary (Ga^ ^In^As P^) alloys. These 
alloys offer an opportunity of engineering several material 
parameters, the magnitude of the fundamental gap being one of 
them. It is possible to tune the band gap of an alloy 
semiconductor by varying the concentration of its components. As 
shown in figure 1.1 the band gap of AlAs is 2.2 eV and that of 
GaAs is about 1.5 eV. By varying x, the concentration of A1 in 
Ga^_^Al^As one can in principle tune the band gap from 1.5 eV to 
2.2 eV. Moreover, unlike amorphous semiconductors (e.g. a-Si:H), 
where it is very difficult to control and regulate the content of 
hydrogen, in these alloys the concentration and the distribution 
of the alloying species can be precisely controlled by using 
sophisticated growth techniques. 

These materials find widespread applications in various 
optoelectronic devices. GaAs^_^Pj^ is used in light emitting 
diodes, Ga A1 As in heterojunction lasers and Hg^ Cd Te in 
infra red detectors. Further the improvement in growth techniques 
such as MBE (molecular beam epitaxy ) or MOCVD ( metallorganic 
chemical vapour deposition) during the last decade has made it 
possible to grow these alloys into super lattices and quantum 
structures. Such structures are important both technologically and 
from the point of view of physics. These structures consist of 
alternating thin layers of different semiconducting materials. The 
material with large band gap acts as a quantum barrier. The height 
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of the barrier can be monitored with the help of alloy 
semiconductors. Thus, these structures lead one to probe the 
quantum regime experimentally. 

The main problem which is encountered in growing high quality 
heterojunctions and super lattices is that the materials of the 
various layers should be properly lattice matched to avoid strain 
in the structure. This severely constrains the choice as can be 
seen from figure 1.1 . Except for the pairs GaAs and AlAs, and 
GaSb and AlSb, for all other pairs of compounds the lattice 
parameter is very poorly matched. Hence the importance of alloys 
like Ga A1 As and Ga^ A1 Sb which are free from the strains of 
lattice mismatch, cannot be overemphasized. 

The primary thrust of this thesis is the study of point 
defects in III - V alloy semiconductors. The alloy semiconductor 
may be regarded as a substitutionally disordered system to a good 
approximation. For example, Ga^ 25 ^® posseses an underlying 
Zinc blende structure. The anion sites are occupied by As atoms 
while each four-fold coordinated As atom is surrounded on an 
average by one A1 and three Ga atoms. One can obtain the 
electronic structure by employing a suitable mean field theory for 
disordered systems. Some of these theories are the Virtual Crystal 
Approximation (VGA) and the Coherent Potential Approximation (CPA) 
In section 1.2 we explain the central idea behind the CPA with two 
simple classical examples. 

The presence of point defects in these alloy semiconductors 
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give rise to very interesting and puzzling physical features, 
metastability and negative U behaviour are some of these. DX 
center and some other defects have been discussed in detail in 
subsequent chapters. Further, unlike amorphous semiconductors, 
these alloy semiconductors provide a possibility of a controlled 
study of defect in a disordered host. 

From the technological point of view the presence of defects 
play a crucial role. They are responsible for device degradation 
(e.g. the dark line defect in the light emitting diodes made of 
Ga^ As P). Point defects in semiconductors also give rise to a 
deep trap inside the forbidden band gap. These deep traps are 
important because they act as recombination centers and control 
the lifetimes of carriers in luminescent devices. They are used to 
increase the probability of radiative electron hole recombination 
in light emitting diodes in order to increase the intensity of 
emitted light. This necessitates that the defect which causes the 
deep trap must be identified so that measures could be taken to 
control its concentration in the semiconductor host. 

The defect level due to an impurity has been obtained via 
the Roster-Slater equation in our calculations. In section 1.3 of 
this chapter we illustrate this formalism by considering a simple 
example. Section 1.4 contains a preview of the various chapters of 


this thesis. 
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1.2 The Coherent Potential Approxination (CPA) 

In subsequent chapters we have used the sophisticated 

Coherent Potential Approximation (CPA) for the description of the 

disordered (alloy) host. The CPA was introduced in condensed 

matter physics exactly two decades ago by Paul Soven^ and David 
2 

Taylor . Its genesis can be traced back to the last century . A 
similar methodology was employed by Bruggeman in the 1930 's to 
understand the conductivity of a random mixture of two metals. 
Seven proposed it as a scheme to understand the electronic 
structure of random binary metal alloys, while Taylor carried out 
a similar program for phonons. Since then the random systems have 
been widely studied and the CPA philosophy has been applied 
extensively to study all types of elementary excitations e.g. 

, 3 

excitons , magnons etc. 

The objective of the CPA is to provide a description of the 

average properties of a random mixture. The CPA formalism can be 

well illustrated by considering two simple classical problems”*. 

Let us consider a linear system of 'N' springs, springs 

of stiffness each and 'N^' springs of stiffness each (N = 

N^+N^) . Let us assume that such a random system can be truly 

represented by another linear system of N springs with, stiffness 

K „„ for each spring, as shown in figure 1.2a. The objective of 
err 

any mean field theory is to express in terms of stiffness of 

two types of springs K , K, and their concentration N /N and N /H 

12 X 

respectively. One choice for 


could be 
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Kgff = (N^/N ) + (N^/H ) (1.1) 

This is known as Virtual Crystal Approximation (VGA) and can 
be compared to Vegard's law (where the bond length of an alloy is 
set equal to the sum of individual bond lengths weighted by the 
concentration of each of its components). Expression (1.1) holds 
only in restricted cases. Here it holds when approaches K^. 

Let us suppose that each of the effective springs stretch by 
an amount x when a tensile force 'F' is applied to the system. 
Thus 

X = \F\/ (1.2) 

or 

N X = N IF]/ (1.3) 

Here 'N x' is the total extension of the system and / N is 

the stiffness of the entire system. 

Now we perform a CPA gedanken . We replace one of the 
effective springs by the first type of spring (figure 1.2b). Let 
x^ be the extension of this spring when a force F is applied to 
it. Hence 


|F1/ Ki 


(1.4) 
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In this series combination of springs we notice that the tension 
is uniform throughout and so the total extension is the sum of 
individual extensions of each spring when the same force is 
applied to each of them separately. Thus, the total extension in 
this case will be 


= (N - 1) X + x^ (1.5) 

Substituting for x and x^ from equation (1.2) and (1.4) we get 

= (N - l)|Fl/ + IF]/ (1.6) 

We repeat the gedanken with the second type of spring (figure 

1.2c). In this case the total extension will be 

X^ = (N - DlF]/ + jFl/ (1.7) 

If springs of stiffness are a true representation of the 

collection of these springs, then the average of the extensions X^ 
and X^ should be equal to the total extension of the effective 
system . 


N^Xy N + N^X^/ N = Nx (1.8) 

X and X from the above equations and 
2 


Substituting for X 
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simplifying we get the following equation 

1 / Kgjj = (N/N) 1/K^ +(N^/N) 1/K^ (1.9) 

In the light of the above problem one sees that unlike VGA 
(where the effective value of a physical parameter is set equal to 
the arithmetic mean of its components ), in CPA one first seeks a 
physical parameter which is additive. In the above problem, 
extension by a given force is such a parameter. The effective 
value of this new parameter (extension) is set equal to the 
arithmetic mean of its components (equation (1.8)). Knowing the 
functional dependence of the earlier parameter (spring constant in 
this case) on this new parameter (extension) one gets the required 
relation (equation (1.9)). Note that for the parallel combination 
of springs, extension for each spring is same and the total force 
applied for a given extension (say unity) is the sum of the forces 
applied to each spring separately to obtain the same (unit) 
extension. Thus, the force per unit extension which is nothing but 
the spring constant, is itself an additive parameter. Hence for 
this arrangement both CPA and VGA give identical results. 

In the above gedanken the substitution of any number of sites 
at any location will yield identical results. It is, however, not 
so easy to find such a local additive parameter for every problem. 

For problems relating to binary alloys, a similar procedure 
of replacing them with an effectively ordered structure is 



1-9 


employed. In the single site CPA one site of this effectively 
ordered structure is removed and replaced by one type of atom. 
This replacement is treated as an impurity in an ordered medium. 
The scattering matrix T is evaluated for this new system. In this 
case the scattering matrix T is additive. Hence the average of 
this matrix over all possible types of insertions is equated to 
zero. That yields the CPA equation. It is, however, observed that 
in this case the results obtained depend on the number of sites 
replaced. The transfer matrix T depends on the environment and it 
is not a local property of a particular atom at a particular site. 
This illustrates the approximate nature of the theory. 

The application of the CPA formalism is not constrained to 
treating alloys and mixtures alone, it can also be applied to 
yield good approximations to many other physical problems. 

Consider the case of Ising spins. The Hamiltonian for this 
system is given by the equation 


n . n 

ti = - J Z a o- 
i- ^ j ^ 

J > 0 

Let 'm' be the average magnetisation per site for this system. In 
the Braggs William Approximation the magnetisation 'm' per site is 
given by the equation 
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m = tanh ftJzm 

where z is the number of nearest neighbours, J is the coupling 
constant and ft iS' the inverse temperature. 

The same expression can be derived in a much simpler way 
using the concepts of CPA. At a particular temperature some of the 

spins point up and others point down randomly (figure 1.3a). We 

substitute this random system by an effective ordered structure 
where the magnetisation at each site is 'm'(figure 1.3b). Now we 
focus on a particular site. We replace this site first by an up 

spin and then by a down spin (figure 1.3c). The change in energy 

of the two system in these two cases will be 

2 

AE = Jzm - Jzm 

up 

AE^ = Jzm^ + Jzm 

down 

But a spin at a particular site can point up or down with a 
Boltzmann probability. Hence the thermal average of the change in 
energy should be zero. So 


-exp((^Jzm( 1-m) )Jzm( 1-m) + exp( -f3Jzm( 1+m) ) Jzm( l+m) _ g 
exp(/3Jzm( 1-m) ) + exp( -rJJzm( 1+m) ) 

The denominator does not blow up because the magnetisation m is 
finite and we are seeking solutions at non zero temperatures. 
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Hence the numerator can be equated to zero, which yields 

(l+m)/(l-m) = exp(2/?Jzin) 
or 

m = tanh ^Jzm 

which is the same as the Braggs William Approximation. 

1.3: Koster Slater Equation 

We have mentioned in section 1.1 of this chapter that the 
deep traps in semiconductors are identified with the help of 
Koster-Slater equation . This is a determinantal equation; 

I I 1 - V I I = 0 

Here is the host Green's function operator. The impurity 

n 

potential operator is 'V'. E, serves to locate the defect level. 

Q 

In general the operators V and are matrices. In this 

introductory presentation we take them as scalars. 

Green's function can be expressed in terms of the host 

density of states' feis follows 


G^(z) 


r n(z ' )d z ' 

J 
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where n(z') is the density of states of the host medium. 

The methodology of obtaining can be well illustrated by 
considering the case of a hypothetical semiconductor whose density 
of states is a sum of two delta functions at energy r and -r. 
Thus , 

n( 2 ) = <5<z + r) + 6(z - r) 

The delta functions at r and -r can be identified with the 
conduction and the valence bands respectively. Figure 1.4 depicts 
sush a semicircular DOS. 

To obtain the energy level due to an impurity of potential 
strength 'V' in such a semiconductor we proceed as follows. We 
first evaluate the host Green's function. Thus, in this case 


G^(2) 


r ^(z'+rid z' 6(z'-r)d z' 
J (z-z') J (z-z') 


= i/(z+r) + i/(z-r) 


Now employing the Koster-Slater equation we get 


2z/(z^-r^) = 1/V 


The above equation is a quadratic function in z, the solution of 
which gives the location of the defect level. Thus, 
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2^ - 22V - = 0 

= z = V i cv^ + 

When V << r 

2 = V - r fi + 

As V approaches zero the defect level tends to - r, as is clear 
from the above equation. Thus as V decreases the defect level 
shifts towards the band edge, for V zero it coincides with the 
bands. When V >> r 


4 * 2 2 ±y2. 

2 = V - V (1 + rVv ) 

2 

From the above equation it is clear that as V increases the term 
2 2 

r /V under the square root sign becomes negligible. Thus as V 
tends to infinity, from either the positive or the negative side, 
one of the solutions for z goes to 2 ero, that is the level shifts 
towards the mid gap as |V| is increased. This phenomenon is known 
as the pinning of the leveland is illustrated in figure 1.4. There 
is, however, a second solution for z in this case which is outside 
the gap and is not of much relevance. 

1.4: Summary 

The examples presented in the previous sections are intended 
to delineate the broad conceptual framework of this thesis. In the 
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following paragraphs we briefly summarize the work presented in 
the subsequent chapters. 

In chapter 2 we present defect energetics calculations for 
some model disordered hosts. We have chosen two types of model 
hosts. One of them is the nearest neighbour one dimensional tight 
binding host and the other is a host that has a semicircular 
density of states. We have studied disorder in these model systems 
using CPA. We have also employed the weak scattering limit 
Nordheim's Correction (NC) for small disorder. The Strong 
Scattering Limit (SSL) has also been used for large values of the 
disorder parameter. The total energy changes in these systems due 
to a point defect has been evaluated using Hellmann Feynman 
theorem. The general trend of the total energy change with the 
defect potential is found to be the same in all the three 
approximations. However the saturation value of the energy change 
for large positive defect potential is underestimated in NC as 
well as in SSL as compared to CPA**. 

Chapter 3 contains an approximate form of CPA developed for 
realistic multiband semiconductors. One of the main results of 
this chapter is that both the CPA and VCA predict identical values 
for a vacancy related trap level. 

In chapter 4 we have concentrated on a particular defect 
namely the DX center which is present in Ga^_^Al^As. This is 
basically a donor related level. However since it is deep inside 
the gap it is believed to be a donor complex and not a simple 
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donor. But there is controversy regarding this assignment. The 
shallow to deep transition of this defect level when the A1 
concentration is increased beyond 20% seems to suggest that the 
defect could be a simple donor, and the disordered environment is 
responsible for its being deep. This idea has been explored in 
Chapter 4 by using the Coherent Potential Approximation (CPA) to 
account for the disorder self consistently. Our results agree with 
the equilibrium Hall measurement^. 

In Chapter 5 we have developed a formalism for helium like 
impurities (the defect atom has two valence electron more than the 
host atom) in th^,, semiconductor matrix. The chemical shift of the 
impurity atom has been taken into account by considering the 
impurity to be inside a cavity of radius r . The influence of the 

O 

impurity within this cavity is unscreened by the medium. We have 
applied this theory to study several classes of defects e.g. (i) 
chalcogen impurities in silicon (ii) anti-site defect in III-V 
semiconductors (a popular model for EL2 in GaAs) (iii) DX center 
(which is viewed by some authors as a negatively charged donor due 
to its EPR invisibility ). 

Finally, in Chapter 6 we have focussed on the analysis of 
experiments. In particular we have concentrated on Deep Level 
Transient Spectroscopy (DLTS) and Paired Temperature Spectroscopy 
(PATS) which are employed to find the location of the gap states 
in semiconductors. We point out that the analysis of the 
results obtained from DLTS and PATS experiments is based on an 
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expression which is valid for the ordered semiconductor. For 
disordered materials the same expression may no longer be valid, 

P lO 

leading to a systematic discrepancy ' . In chapter 6 we have also 
focussed on the effect of disorder on M center, a metastable 
defect in InP and its alloys^^. 
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FIGURE CAPTIONS 


Figure 1.1: A plot of band gap vs. the lattice parameter for some 
typical III-V compound semiconductors. 

Figure 1.2: The gedanken (thought experiment) described in section 
1.2, (a) depicts a linear array of effective springs, (b) and (c) 
represent the act of removing one of the effective springs and 
replacing it by the springs of the first and the second kind 
respectively . 

r r f' 

Figure 1.3: The gedanken on Ising system described in section 1.2, 
(a) shows the random arrangement of spins at temperature 'T', (b) 
shows the effective system which replaces the random arrangement 
of the spins, (c) depicts the act of removing one of the effective 
spins and replacing it by an up spin and a down spin respectively. 

Figure 1.4: An idealized semiconductor DOS. Both the valence and 
the conduction band are represented by a <5 function of equal 
strength. The resulting defect level is shown below. 
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CHAPTER 2 : DEFECT EHERGETICS : MODEL ANALYSIS 

2.1: Prelininary remarks 

In this chapter we present an analysis and our results of 
total energy calculation using Hellmann Feynman theorem for a 
substitutional defect in a model disordered system. This analysis 
could be useful for the study of those alloys which are 
subst itut ionally disordered. A detailed description of the 
treatment of a realistic alloy semiconductor will be given in 
subsequent chapters. 

The description of the host medium which is disordered, 
entails the use of mean field theory. There are several mean field 
theories. A brief recapitulation of these mean field theories 
could be useful in this context. 

Consider a single band, orthogonal, tight binding, 
Hamiltonian with diagonal disorder due to the presence of two 
types of atoms A and B. 

n . n 

H = Z £ |m> <m| - t 51 |m> <n I (2.1) 

m n 

Here ^ ^ with a probability C^ 

m A .A 

= £ with a probability 

Here C and C are the concentrations of the two types of atoms. 

A B 
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|in> and |n> denote the Wannier states at the sites m and n 
respectively. The transfer term connects the nearest neighbours 
only . 

In the mean field approximation c is replaced by a 

m 

renormalized self energy <^(z) , where z is the electron energy- H 

D 

then assumes an effective ordered form. 

n . n 

m 1®^ - t 2 lm> <n| (2.2) 

m n 

An ordered Hamiltonian on the other hand has the form 

n . n 

~ ^ m 1®^ '^®l " t 2 |m> <nl (2.3) 

m n 

Here u is independent of the electron energy z . Since is just a 
constant in this case, the zero of the energy can always be set 
such that jJ is zero. Hence 

H = H + (o'(z) - ^ Imxm] 

or H = H + <3'(z) 2 Inixmj (2.4) 

The diagonal part of the Green's function at the site labelled 
zero is 


F (z) = <0 1 / (z-H^) |0> = 1/N 2: (l/(z-£(k)) 

o o , 


(2.5) 
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F( 2 ) = <01 l/(z-H) 10> = 1/N E (l/(z-e(k)-<y(z)) (2.6) 

k 

Here ^(k) is the energy eigenvalue of the ordered system. It is 
clear from equations (2.5) and (2.6) that the diagonal part of the 
ordered and disordered Green's function are related by the 
following equation 


F(z) = Y^iz-cr(z)) (2.7) 

The self energy oiz) is given in the various mean field 
approximations as follows: 

(i) The Virtual Crystal Approximation (VGA); In the VGA <^(z) is 
defined as the simple average 

<y( z") = G £ + G £ - £ (2.8) 

A A B B 

(ii) The Coherent Potential Approximation (CPA): In the more 
sophisticated CPA it is obtained by the following self consistent 
equation 


<y(z) = -e - (£ - cf(z')') Ffz) ( £„ - c/iz)) 

A ® 


(2.9) 


There are two limiting cases of CPA which are of interest to us 
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(a) The strong scattering limit (SSL) where 

c'(z) = £ + C C /(7.-£' ) (2.10) 

A B 

Here s' c + C £: 

A B BA 

& u ^ - s )/ 2t 

A B 

In the SSL u is presumed to he ler^fe 

(b) The weak scattering limit, Nordheim's Correction (NC)^ 

o-(2) = ^ + (2.11) 

Here ~ 

The density of states (DOS) is related to the imaginary part 
of the Green's function. Hence the DOS of the system can be known 
from the knowledge of the Green s function. 

If a single site impurity is added to a medium the potential 
term in the Hamiltonian changes. The impurity potential is defined 
as the difference between the total Hamiltonian with impurity and 
the host Hamiltonian. 


V = H - H, (2.12) 

Imp r 

This impurity potential perturbs the host electronic structure. It 
may give rise to a defect level inside the band gap. The 
conditions for the presence of the defect level inside the band 
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and its location will be discussed in section. 2.2. The 
perturbation of the host DOS and the presence of a defect level in 
the gap results in changes in total electronic energy. In section 
2.3 we discuss a formalism to calculate the total energy changes 
using Hellmann Feynman theorem. Section 2.4 contains the results 
of the calculation of total energy changes in two types of model 
disordered systems. 

2.2: Determination of the Defect Level in the Gap. 

Consider the ordered, one dimensional, tight binding, 
Hamiltonian of equation (2.3) in which one atom of different 
species replaces one of the parent atoms. This atom is then 
referred to as a substitutional impurity. The Hamiltonian of the 
system can be written using equation (2.12) as 


Here V is the impurity potential and H, is the host 

imp • ^ 

Hamiltonian. Assuming that the substitutional impurity is neutral 
and unrelaxed, we obtain 

V. = |1> (^ - ^^)<l! = 11> V <1| (2.13) 

Here '1' is the defect site. Since we have assumed the impurity to 
be unrelaxed, that is, there is no rearrangement of neighbouring 
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atoms. Hence 

there is no 

change in 

the 

off-diagonal 

term of 

the 

Hamiltonian . 

The Green's 

function 

of 

this system 

is given 

by 


(Economou 1979)^ 


or 

where 


G(z) = 1/(2 - H) 


= 1/(2 - H. - V. ) 

h tmp 

G(2) = G^(2) + G^(Z) V, G(2) 

h h tmp 


G^(2) 


l/(z-H^) 


G.(z) is the host crystal Green's function. Defining the 

n 

scattering matrix T as 


T 


V. / 

imp 


(1 - G^(2) 


V. ) 

tmp 


It can be shown that 


G(z) = G^(2) + G^(z) T G^(z) 


T can also be written in an expanded form 


T 


V. + V. G 

tmp tmp 

Vllxli + 
11>V<1 I / Cl 


1 1 ><1 iG^(z) 1 ixl 1 

- G^d.l.z)) 

n 


+ 


where 
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G^(l,l,z) = <l|G^(z)ll> 

Hence 

G(2) = G^(2) + G^(2)ll>V<llG^(2)/(l - G^(1,1,2)) 

The poles of the Green's function give the energy levels. From the 
above equation we see that the impurity Green's function is the 
sum of two terms. The first term is the host Green's function and 
the second term is a function of the impurity potential. The poles 
of the second term give additional levels over and above the host 
crystal levels, at say E which is given by the solution of the 

p 

following equation. 


1 - V G.(1,1,E^) = 0 

h p 

or 

G^d.l.E ) = 1/V (2. 14) 

h p 

3 

The above equation is known as the Roster Slater equation for the 
single band case. It is, however, not necessary that E^ should 
always be outside the band. We may have a host DOS such that for a 
small value of the impurity potential there is no solution of 
equation (2.14) outside the band. This situation is labelled as 
resonant. Inside the band the Green s function is complex and the 
T matrix could be written as 

T = V 1 1x1 1 / (1 - V G^(l,l,Ep) - iV Gj^(l,l,Ep)) 
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The superscripts R and I indicate the real and imaginary part of 
the Green's function. Further, 


= 1/V 

yields the resonance state inside the band gap. 

It is observed that for some host DOS there is a threshold 
value of impurity potential. For values of defect potential larger 
than the threshold a level appears in the gap. It is easy to see 
from the K-S equation that for the threshold to be zero the real 
Green's function just outside the band edge should be infinite. 

2.3: Determination of Total Energy Changes 

The gap level by itself does not present the complete picture 
of a defect in a semiconductor. The total change in energy is 
evaluated not merely by considering the gap level, but also the 
perturbation of the host density of states and the location of the 
fermi level (z^.). This total energy is related to the heat of 
solution . 

In this section we present the formalism to employ the 
Hellmann Feynman theorem to obtain the total energy change when 
the defect is introduced into a model disordered host. 

The Hellmann Feynman'^ theorem states that if H(X) be a 
hermitian operator', which depends on the real parameter. ' and 
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|V'(>^)> a normalized eigenvector with eigen-value E(X) then 


dX 


(2.15) 


The total energy of the system is given by 


Ef 

E = J w p(u>) dco (2.16) 

-00 


where p(<*i) ; density of states 

: # of states per unit energy 


But p('^) can be expressed in terms of the Green's function 


= 1 /n Z Tm (T-, , (w) 


Here *3-j^j^(co) is the diagonal part of the Green's function 
The total energy changes can be expressed as 


Ef 



r 

C f 'j ' 

t 


AE = - ’ 

^ l/n 1 00 Im ■< 

lGimp(05)J^^ - G^^(OO) 

> doo 

(2.17) 


-00 


Here (G ) is the diagonal part of the Green's function of the 

Imp ii 
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host with the impurity atom and G^^(co) is that of the host medium. 
If the host medium is ordered then G..(w) is same for all 'i's 
because of translational invariance. For disordered host the 
medium is not translationally invariant, however, as mentioned in 
section 2.1 the medium can be replaced by an effective ordered 
system. Thus for this case too the diagonal host Green's function 
will have the same value for all sites. 

Since the impurity atom is located at a particular site, 
(G. ).. is site dependent. Thus one needs to evaluate (G. ).. 

vmp <.». ' vmp V.V 

for every site separately, which is not easy. This difficulty is 
overcome by taking recourse to the Hellmann Feynman theorem. The 
Hamiltonian of the host system is 


H 


h 


a Z 

m 


n . n 

|m> <mj + t Z jm> <nj 
m n 


and that with an impurity is 


•m 

n . n 

H = O’ Z: lm> <m I + t ^ jm> <n | + (£ - o') 11><1| 

m ^ n 

n . n 

=: <y H |in> <in| + t ^ |ni> <n| + V |1><1| 

m ’ ‘ ^ ‘ 


Here the Hamiltonian depends on the parameter V, the impurity 
potential. Using Hellmann Feynman theorem one gets the following 
equation as shown by Villasenor et al 
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Ef 

J <V'|1> <ll\y> dw 
00 
Ef 

J 1 <^1 1> dw 
00 
Ef 

J P/,w,V) dw (2.18) 

00 

P^(<^»V) is the local density of states at the impurity site and is 
given in terms of the diagonal part of the impurity Green's 
function as follows 

p^(<o,V) = -1/n Im (Gtr«p(<o,V))^^ 

The evaluation of the diagonal part of the impurity Green's 
function at the defect site is straightforward. It is given by the 
expression 

Integrating equation (2.18) with respect to the impurity potentia 
the total energy change can be easily obtained. 

2.4: Results for Total Energy Change in Two Model Cases 

We investigate the effect of various mean field 
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approximations on the total energy change for two types of model 
disordered hosts. 

(i) The first model is that of one dimensional tight binding 
system . In this case the diagonal part of the Green's function is 
given by 

F(2) = (2.19) 

(ii) The second model is that of a realistic three dimensional 
semicircular density of states where 

F<z) = 2 ((z-o'Cz)) - (z-o-( 2 ))^ - (2.20) 

Our calculations are for a 50-50 alloy. We will focus on two 
different physical situations, one in which the atoms of the alloy 
are chemically very different from each other, and the other in 
which they are almost similar. These situations are characterized 
by large and small values of the disorder parameter 'u'. The 
disorder parameter 'u' is defined as the difference between the 
diagonal energies of the two atoms divided by the off diagonal 
term ' t ' . 

u = -e ) / 2t 

A o 
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We next present results for the change in total electronic 
energy due to a single site impurity in the above mentioned model 
disordered hosts. 

In Figure 2.1 we predict the change in electronic energy 

as a function of the defect potential for various band fillings 

for the realistic semicircular host . The Fermi level is denoted 

by Zj. Note that z^ = 0 corresponds to the half filled band, z^ = 

2.5 (-2.5) correspond to approximately three quarter (one quarter) 

filled bands . Note further that V - £ . . The disordered host is 

imp 

described in the Coherent Potential Approximation . The disorder 
parameter u is large ( = 5.0) . The defect potential (V) perturbs 
the host band as well as creates a localized level outside the 
host band. For negative values of the defect potential the 
localized states appear below the band, consequently lowering the 
energy . From the figure it is clear that for large negative 
defect potentials the plot is linear irrespective of the band 
filling . The position of the defect level is proportional to V 
for large V .This being the dominant contribution to the energy , 
accounts for the linearity of the plot at large defect potentials. 
The level characteristic are the same for positive V , this time 
however the level being above the band, it is left unoccupied. 
Consequently we expect AE to approach a saturation value as V is 
increased . 

Figure 2.2 represents the same plot as that of figure 2.1 
with identical parameters, except that now the disordered host is 
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described by employing the Nordheim's correction . The general 
behaviour of the curve is similar to figure 2.1 (e.g linearity at 
large negative values of V and saturation at large positive V). 
The saturation occurs at approximately the same value of V as that 
in figure 2.1 . However the saturation value of AE is considerably 
lowered ( approximately 15% ) as compared with the results in 
figure 2.1. 

In Figure 2.3 we have used the strong scattering limit 
to describe the host . The general trend of the curve is similar 
to that of CPA . However, the saturation value for large positive 
V is once again underestimated as in the Nordheim's correction. 

A similar analysis was repeated for the same model host 
which was weakly disordered (e.g. u = 0.8). The results obtained 
by employing Nordheim's Correction were found to be in very good 
agreement with the CPA. 

In figure 2.4 and 2.5 we investigate the effect of two 
mean field theories: CPA and NC respectively on the defect energy 
in a one dimensional host .The value of the disorder parameter is 
chosen to be the same as in figure 2.1 and 2.2 ( u = 5.0) The 
overall behaviour is similar to figures 2.1 - 2.3 (e.g. linearity 
at large negative values of V and saturation of AE at large 
positive values of V). We notice however that the Nordheim s 
Correction severely underestimates the saturation energy AE for 
large positive V .In the strong scattering limit the results seem 
to improve, but nevertheless saturation value differ considerably 



2-15 


from the Coherent Potential Approximation 

For weakly disordered one dimensional host ( u = 0.8) 

the results obtained by using Nordheim's Correction was found to 
agree very closely with the Coherent Potential description. 

The results of this chapter have been presented in detail in 
Reference 6. 
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FIGURE CAPTIOHS 

Figure 2.1: Electronic energy change AE vs. the defect potential V 
for various band fillings . Note that = 0 corresponds to a half 
filled band . The model considered is the realistic three 
dimensional host . 

Figure 2.2: Parameters and the model is the same as figure 2.1 
except that the Nordheim's Correction is used to describe the 
host. Note the general similarity of the figures 2.1 and 2.2. 

Figure 2.3: Electronic energy change AE vs. the defect potential V 
for various band fillings . Parameters and the model used is same 
as figure 2.1 except that the Strong Scattering Limit is used to 
describe the host .Note that there is no improvement over the 
Nordheim's result. 

Figure 2.4: Plot of electronic energy vs. the defect potential for 
various band fillings for a one dimensional host . The host is 
described by the Coherent Potential Approximation. 


Figure 2.5: Same as figure 2.4 except that the host is described 
by the Nordheim's correction . Note that NC proves to be a very 
poor approximation. 
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CHAPTER 3 ; MULTIBAND CPA FOR REALISTIC 
SEMICONDUCTOR ALLOYS 


3.1: Preliminary Remarks 

In this chapter we will study defect energetics in realistic 
semiconductor alloys. We divide the framework necessary to study 
electronic energetics associated with defects into various stages. 
The first stage is the determination of the host electronic 
structure. For this an alloy semiconductor may be regarded as a 
substitutionally disordered system. For example 

posseses an underlying zinc-blende structure. Each four-fold 
coordinated As atom is surrounded on an average by one A1 and 
three Ga atoms. Since the cation site can be occupied by either A1 
or Ga with probabilities which depend on their respective 
concentrations, the analysis of chapter 2 can be employed. This 
disordered structure, in which the occupancy of the cation site is 
probablistic, can be replaced using MET by an effective ordered 
structure in which all cation sites are considered to be 
equivalent. This replacement is, however, at the cost of several 
physical characteristics of the original system. In the original 
disordered system there can be localization effects, which are 

u A ¥ 1-hp pffective ordered crystal which replaces 

largely washed out by the errecci 

it . 

4 - * ^ disordered medium could be 

A realistic description of a aiooru 

1 rto oliT^iter calculations averaged over all 
achieved by very large clus 
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possible configurations. This method is computationally feasible 
only for very small cluster sizes, which probably do not do 
justice to real systems. In our work we have concentrated on the 
mean field description of the disordered medium. 

The second stage is to determine the mean defect level in the 
gap Z^(V). Here V is the defect potential. The mean gap level 
Z^(V) is then given by the well known Koster Slater equation^. 

In section 3.2 we present the LCAO tight binding scheme to 
calculate the electronic structure of the host. In section 3.3 we 
introduce the Coherent Potential Approximation (CPA) formalism and 
describe the methodology employed to incorporate the CPA for the 
multiband case. Section 3.4 contains the description of a proposed 
deep level approximation. Section 3.5 contains results and 
discussion . 

3.2: LCAO Tight Binding Scheme 

As noted in the earlier chapter, the host electronic 
structure should be known in order to predict the presence of gap 
level due to an impurity. In this section we describe the method 
adopted to find the host structure. 

We employ the empirical tight binding scheme for the 
evaluation of the host band structure. We wi’ite the crystal wave 
function as a linear combination of Bloch functions |nbk> each of 
which depends on only one type of orbital. Thus our basis set is 


given by 
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lnbk> = N Z exp( ikR. + ikx )lnbR.> (3.1) 

V ,b ' b ' L 

|nbR^> are the atomic orbitals centered on the N crystal sites. 
The wavevectors k lie in the first Brillouin zone of the 
reciprocal lattice. The band index is n and b is the sublattice 
index. The primitive vectors of the (fee) lattice are 

= a/2 (x + y) 

Aj, = a/2 (y + z) 

A^ = a/2 (x + z) 

Here 'a' is the lattice constant and x, y, z are unit vectors in 
the Cartesian system. The primitive vectors of the bcc reciprocal 
lattice are 


a^ = In/z. (x + y - z) 

a^ = 2n/a (z + y - x) 

a^ = 2n/a (x + y + z) 

For a bcc reciprocal space the k vectors are defined by 

k = 2n/a 

Here each of the M ,M, and M. is less than or equal to unity. 
Furthermore, they satisfy the condition 

M^+Mz+M^ ^ 3/2 

The above two conditions on M^'s define the unit cell (the first 
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Brilloin zone) in k space (bcc type). 

The zinc-blende structure comprises of two interpenetrating 
fee lattices one for the cation and other for the anion. The index 
b in equation (3.1) denotes whether the site is an anion site 
or a cation site. Assuming that the anion positions are given by 
R 's . We define 

I 


to take the cation site into account. 

The valence electrons are the ones which participate in the 
covalent bonding. The compounds with which we will be dealing in 
our calculations have zinc-blende structure and their valence 
orbitals comprise of one s and three p shells. The index n in 
equation (3.1) stands for the s, px, py, and pz orbitals. In the 
basis |nbk> the Hamiltonian matrix is of size (8x8). The matrix 
elements of the Hamiltonian are obtained empirically by matching 
them to the experimentally obtained band structure at high 
symmetry points. However, attempts to fit the conduction bands 
fail. To match the conduction bands two excited s orbitals, one at 
the anion and the other at the cation site are introduced in the 
basis set of ]nbk>, thus rendering the Hamiltonian matrix to be of 
size (10x10). 

A representative table of Hamiltonian matrix elements is 
shown by Harrison It has been assumed that the Hamiltonian couples 
only the nearest neighbours. The empirical values of the diagonal 
and the off diagonal terms E(njb) and have been 
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3 

listed by Vogl et al for several semiconductor elements and 
compounds . 

In the case of semiconductor alloys, there is a substitional 
disorder either at the anion or cation site (for ternary alloys) 
and at both sites (for quarternary alloys). We have used the VGA 
approximation discussed in chapter 2 to determine the input matrix 
elements of alloys (Ho et al 1983)“*. The CPA corrections to these 
is described later in section 3.3. 

Consider an alloy of composition A B CD - Let us denote 

X dL—x Y 1“Y 

the lattice constant (or bond length) for the compounds AC,AD,BC 
and BD by di, dz, da and d-* respectively. The representative 
diagonal matrix elements of the above four compounds are written 
as (AC)nn, (AD)nn, (BC)r.n and (BD)r.r.; and the off diagonal terms 
as (AC)n n , (AD)r. r, (BC)n r. and (BD)n n 

xZ X Z jLZ X ^ 

The diagonal matrix element of the VGA averaged Hamiltonian 
for the alloy is given as follows 


(H)r>n = xy (AC)nn + x(l-y)(AD)nn + (l-X)y(BC)nn + ( 1-X ) ( 1 -y ) ( BD ) nT^ 

(3.2) 


(H)n n 

12 


= xy(d/d)^(AC)n^n^ + x(l-y)(d 2 /d)^(AD)n^n^ 

+ (l-x)y(d/d)^BC)n^n^ + ( l-x)( 1-y )(dyd )^BD )n^n^ 


(3.3a) 


where 

d = xyd^ + x(l-y)d 2 + (l-x)ydg + ( l-x)( 1-y )d_^ 


(3.3b) 


Notice that the above equation respects the inverse square scaling 
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with the bond length of the off diagonal element of the 
Hamiltonian as enunciated by Harrison( 1980 ) . 

D iagona 1 is ing the Hamiltonian yields the band structure of 
the system. We, however, need to find the Greenes function of the 
system. The Green s function can be evaluated by the knowledge of 
the density of states. The calculation of the density of states 
from the band structure was premised on the methodology proposed 
by Lehman and Taut'* in which the Brilloin zone is divided into 
tetrahedrons. A detailed description of this scheme is given by 
Rath and Freeman( 1975)**. 


3.3: Multiband CPA Formalisn 

As mentioned in the earlier chapter the CPA is a better 
method of treating disorder. We will illustrate in this section 
the method adopted by us to introduce a CPA description of the 
disordered host. 

We are dealing with substitutional and not positional 
disorder of the lattice. Hence, the underlying symmetry of the 
lattice remains the same as that of the tetrahedral group (T^). 
The basis functions which we have chosen for our calculations 
decompose into four irreducible representations of the T^ group, 
two are of A^ type (of dimensionality one each ) and the remaining 
two are of the type T^^ (dimensionality three each). A^ is a 
totally symmetric representation which is invariant under all 
group transformations. We will refer to it as an S state. The 
irreducible representation T^ has the same transformation 
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properties as the simple functions x, y, and 2 under the T, 
symmetry operations. The basis functions for this irreducible 
representation could be chosen a.s if> , y/ , and which have an 

y z 

axial symmetry about OX, OY, and OZ cartesian axes respectively. 
We shall refer to them as P states. 

The Hamiltonian expressed in the above symmetry adapted basis 
will be seperated into two blocks as follows. 




H 


Ha^ 0 
0 Ht 

2 


( 3 . 4 ) 


Ha^ is a matrix of size (2x2) corresponding to two 

representations. corresponds to the two T^ representations and 

is of size 6x6. The two representations of the A as well as the 

1 

T^ have, as their basis, functions which are centered either at 
the anion or at the cation site. The Hamiltonian further does not 
couple the basis functions of T^, centered at the same site but 
having axial symmetry about a different axis. Thus Ht^ reduces to 
three (2x2) matrices corresponding to basis functions having 
symmetry about x, y and z axes respectively. These three sets of 
matrices however are identical since the tetrahedral structure 
does not distinguish between the three axial directions. 

A disordered Hamiltonian can be broken into a random part 'V' 
and a configuration independent part W. Thus 


H = W + V 


( 3 . 5 ) 
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The random part is assumed to be site diagonal hence 

H=W + 2Vj (3.6) 

j 

Here 'j' represents the site index. As mentioned in section 3.2 we 
have chosen W to be our input VGA Hamiltonian, and we will be 
incorporating the CPA corrections due to the random part SVj to 
this. The random potential can be written as 

V = r lmR.> {£. - <mR.i (3.7) 

^ J J ’ 

Here m stands for a S like state and three P like states. ^ 

takes values £ ox e depending on whether the atom A or B is 
occupying the site, s is the VGA energy defined in the earlier 
section for W. 

The CPA formalism replaces the configuration averaged one 
electron Green's function by an effective Green s function. 


G 


«ff 


1/(2 - H ,„) 


H ,, is defined as 

off 


H , = w + 

eft 


eff 


= » - 2 b 


(3.8a) 

(3.8b) 


T.. is termed 

j 


as the self energy of the system and 


is obtained 


from 
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the condition such that the single site scattering matrix t when 

j 

averaged over all configurations is zero. E can be viewed as a 
block diagonal matrix, each block is identified by a particular 
site j and is of size (8x8) for the above mentioned eight 
states. Furthermore, the translational symmetry of the averaged 
crystal ensures that each block is identical. Thus Bloch's theorem 
can be invoked and E can be represented in a reduced Bloch space 
of size (8x8). The tetrahedral symmetry further decomposes E into 
two blocks of and T^^ representations. We write E as 


E 


0 

0 Ex 


Ea = 
1 


0 

0 Ea^ 


where the superscripts 'a' and 'c' stand for anion and cation 
respectively , 


Zt 

2 X 


Zt 


2 


Zt 

2y 


Zt 


2z 


Zt . Zt can be further decomposed into cation 

2x ^ 2y 2z 


Each of the Zt 
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and anion subspaces as has been shown for The off diagonal 
elements of 21 have been neglected, which is in spirit with the 
approximation invoked by Chen and Sher^. 

Symmetry about x, y, and z axis further ensures that 


Zt = Et = Zt - Zt 

2x 2y 22 2 

The CPA equation (see equation 2.9) for this multiband case is a 
matrix equation of the form 

Z( 2 ) = -(1^ - Z(z))F(z)(£^ - Z(z)) 

The above mentioned symmetry reduces it to a set of four coupled 
scalar equations 




- 


r^d 




r,a 




r,d 




(3.9) 


Here 'r' refers to the symmetry and 'd' to the cation or anion 
subspaces. These equations are coupled via the Green s function. 
An approximate method of decoupling the above equations can be 

B 

devised using the partitioning technique of Lowdin . 

The component Green's function ^ can be evaluated from the 
component density of states as follows. 


F 


r,d 


= <0 r,d j r,d> 
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Here 'j' is the site index. Since all sites of the effective 
system are identical, we can express F , as trace over all sites 

r,a 


F ^ = l/N Tr <r,d 

r,a 


I r,d> 


N 


I 


<r ,d 


n , k 


( Z 


n,k >1^ 

- > 


[ 


n"(2-.)*^ I I 


l<r,d I n,k>|^ d^k 
^ ^ - ''r,d - > 


_ if d£ 

- N J (2 - ^ ^ ) 


tY l<r,d|n,k>|^<5(£-e (k)) d^k 

N(2n)^ ^ J 


The factor in the square bracket is the component density of 
states p(r,d). Thus 


F 


r,d 


1_ r ^r,d^^> 

N J 2 - - - 


(3.10) 


Using equations (3.9) and (3.10) the CPA corrections to the 
Green's function can be determined. 


3.4: Point Defects in Multiband Systems 

In this section we focus on point defects. These could be 
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either substitutional impurities or vacancies present in the 
disordered alloy semiconductors. The description of the host 
structure has been described in detail in the earlier sections of 
this chapter. 

The defect potential caused by the substitution of an 
impurity atom for a host atom has the following three distinct 
contributions (i) a short range central cell potential (ii) a long 
range Coulomb potential (iii) an electron-lattice interaction 
caused by the lattice relaxation around the impurity. We will 
ignore the screened Coulomb potential and assume a simple 
realistic model in which there is no lattice relaxation. We will 
deal essentially with the central cell potential which is the 
difference between the atomic potentials of the impurity and the 
host atom. 

The single site defect can be identified with of equation 
3.6. In this case the index ' j ' corresponds to a particular site 
and is not a running index over all sites. 

Thus the defect potential can be made block diagonal in this 
symmetry adapted basis; it contains one (2x2) a^ symmetry matrix 
and three (2x2) t^ symmetry matrices. The three (2x2) P type 
symmetry matrices will be identical. The defect potential is then 


V/ 


r 

c 

c , t mp 


e 


r 

c , host 


V. " - v/. 


V. ^ - v," 


r 

e 

a , t mp 


<x , host 
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r represents the symmetry and i & j r {o, a} correspond to 
cation and anion respectively. 

We approximate the atomic energy difference between two atoms 
in a solid to be similar to the difference between free atoms. We 
also neglect the off-diagonal elements as they scale as the 
inverse square of the bond length and we have accounted for these 
in equations (3.2) and (3.3). Thus, if the impurity is at the 
cation site then 


v/. 


r r _ 

^ ^ 0 
c , mp c , hos I 

0 0 


The bound state energy E. can be obtained by the solution of the 

O 

following K-S type equation. 

3.5: Deep Level Approxiuation 

If the defect level is a simple vacancy on a cation or an 
anion sublattice then the corresponding defect potential is 
infinity. From equation- (3.11) we see that the defect level due to 
vacancy will then be given by the zeroes of the Here r 

denotes the symmetry of the level and d the sublattice which has 
vacancy as a defect. 

Substituting the value of in equation (3.9)we get the 



value of 21 
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r,d' which is zero. One thusobserves that for the 
vacancy level the CPA result is identical to the VGA. This 
suggests that a perturbation approach to an ultra deep level in 
alloy semiconductor can be employed. We call this the deep level 
approximation (DLA). This DLA is obtained from the self consistent 
CPA equation by retaining first order terms in inverse defect 
potential ( 1/V) . 

The self consistent CPA equation is (dropping the subscript 
for lattice and symmetry ) 

Using the K-S equation 


F(zd) = 1/V = <5 


where <5 is a small parameter. Thus 


Z(2^) = - Z(Z^)) 6 - ECz^) 


Retaining first order terms in <5 one gets 


E(z^) 


-6 


s 

A 


e 

B 


Moreover 




(3.12) 
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where 

= 2- - 


or 


A B 


(3.13) 


Let 2 ^ be the location of the defect level for the vacancy. 
Expanding ^ about one gets on neglecting higher 
powers 


Fvca<2'> = 




(z 


■2c>^vca'^2c> 


We have chosen <5 to be sufficiently small hence z' can be assumed 
to be quite close to z . But F,, (z ) = 0, since z is by 

C VGA c C 

definition a vacancy level and the VGA result is exact for it. 

Hence 

6 = (z'-z ) F' (z ) (3.14) 

' c' VCA^ ^ 

From equation (3.13) and (3.14) 


= 


^/^'vCA<^c> + ^A 


(3.15) 


Thus knowing the vacancy level and the slope of the VGA Green s 
function at the vacancy level the deep defect level for a given 
defect potential can be directly read from equation (3.15). 



3-16 


3.6: Results 

In this section we present some typical results of our defect 
level calculations . 

Figure 3.1 depicts the t^ cation defect level in the 
fundamental gap of GaAs^ The cation vacancy level for this 

system is at 0.33eV above the valence band edge. The agreement 
between DLA and CPA is excellent near z 

vac 

Figure 3.2 depicts a similar comparison between the DLA and 
the CPA for the anion defect level in the fundamental gap of 
Ga^ results deviate from the CPA as we move away 

from the vacancy level which is at z = 1.35eV. Figure 3.3 shows 

vac 

a t_ cation level in Ga A1 As - 

As demonstrated in section 3.5, the VCA and CPA give 

identical results for the vacancy level in an alloy semiconductor. 

It is worthwhile exploring the differences between the two for a 

general gap level. It must be noted that when the host energies 

for the two species are close to each other, then the VCA ought to 

be satisfactory (small scattering regime). Let us consider the 

case of GaAs Sb • From Vogl^ parameters Table we observe that 
0.8 0.2 

the cation p-orbital energy E(p,c) for GaAs and GaSb differ only 
by 20%. Thus in figure 3.4 , the defect level as charted out in 
the CPA is almost indistinguishable from the VCA. 

When the energy difference between the two species is large 
we expect CPA and VCA to differ considerably. That is borne out by 
a plot of figure 3.5 in which the a^ cation defect is plotted for 
Ga A1 As- The energy E(s,c) differs for two systems GaAs and 

1 — X X 
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AlAs in this case by 56%. 

In section 3.3 we have proposed a simplified multiband CPA 
formalism. It is possible to generalize it. An interesting 
generalisation using the Effective Medium Approximation (EMA) of 
Roth^'*^*^ has recently been proposed. Some of the results of this 
chapter is presented in reference 11. 
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FIGURE CAPTIONS 

Figure 3.1: A Plot of the inverse defect potential strength vs. 

the t cation defect level in the fundamental gap of GaAs Sb 

The solid and the broken lines show the results obtained by using 

CPA and DLA respectively. Note that the vacancy level shown by an 

arrow is predicted at 0.33eV by both the theories. 

Figure 3.2: A Plot of the inverse defect potential strength vs. 

the t, anion defect level in the fundamental gap of Ga Al^ As. 

The solid and the broken lines show the results obtained by using 

CPA and DLA respectively. Note that the vacancy level shown by an 

arrow is predicted at 1.35eV by both the theories. 

Figure 3.3: A Plot of the inverse defect potential strength vs. 

the t^ cation defect level in the fundamental gap of Ga^_ ^Al^^ ^As . 

The solid and the broken lines show the results obtained by using 

CPA and DLA respectively. Note that the vacancy level shown by an 

arrow is predicted at 0.39eV by both the theories. 

Figure 3.4: A plot of defect potential strength vs. the t^ cation 

defect level in GaAs Sb as charted out using CPA( ) and 

0.0 o . z 

VGA ( ). The common vacancy level predicted by the theories is 

shown by an arrow at 0.327eV. 


Figure 3.5: A plot of defect potential strength vs. the a^ 


cation 
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defect level in 0:^0. 4 ^ ^o. charted out using CPA( ) and 

VGA ( ). Note there is no vacancy level within the fundamental 

gap . 
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CHAPTER 4 ; DX CENTER: SUBSTITUTIONAL DONOR DEFECT 
IN ALLOY SEMICONDUCTOR 

4.1: Introduction 

The DX center has attracted the attention of a large number 
of semiconductor physicists on account of its anomalous 
properties. The center, present in the alloy semiconductor 
G’a^_^Al^As presents a challenge both from a basic physics 
and a technological point of view. 

The center as its name suggests is related to a donor 
impurity in Ga^_^Al^As. Group VI impurities (S,Se,Te) as well as 
Group IV impurities (Si,Ge,Sn) in Ga^_^Al^As give rise to this 
defect. The concentration of this defect is approximately equal to 
the donor concentration (N N ) and is independent of the 

DX Donor 

method or condition of epitaxial growth. 

The DX center exhibits divergent activation energies when 

probed thermally and optically, the latter being anomalously 

large . The energies reported depend also on the time scale over 

which the experiments are performed. The Hall activation energy 

is distinct from the Deep Level Transient Spectroscopy (DLTS) 

activation energy E . The former is a slow eQuilibrium measurement 

e 

while the latter is a relatively 'fast non equilibrium process. 
Further, there is a barrier to capture which leads to 

persistent photoconductivity (ppo) and consequent parasitic 
effects in the devices. Lang^ has reviewed some of the salient 
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experimental features of the DX center recently. 

This center has been reported in 

(i) Ga^_/l^As (0.2 < X < 0.45 ) 

(ii) GaAs under pressure 2: 20 kbar . 

(iii) GaAs under heavy doping ( N > 10^^ p«r cm^ ) 

Donor 

Theories proposed to account for its peculiar characteristics 
fall into three categories. 

z 

(i) The large lattice relaxation model of Lang et al . 

(ii) The effective mass theory where the donor is presumed to be 
coupled to an energetically higher symmetry point (L or X) in the 

3 

conduction band . 

(iii) The alloy model where the substitutional disorder in 
Ga Al As is held to be primarily responsible for its anomalous 

1-X X 

behaviour^ . 

It has besn widely conjectured that the DX center is a donor 
vacancy complex. There is evidence however that the DX center is a 

5 

simple substitutional donor and not a complex. Mizuta et al have 
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argued for this model based on Hall effect and DLTS experiments. 
Morgan has proposed that it is a simple substitutional donor 
which may be displaced albeit slightly from its normal centered 
position. Arguments based on coupling to L, X and r points all 
presuppose a substitutional impurity. 

Our investigation is limited to this 'substitutional* 
premise. Within this paradigm we shall present calculations which 
are non perturbative and which self consistently account for the 
alloy disorder. The formalism is based on the tight binding 
Green's function method. It is well suited for cases where the 
defect is localized. 

In section 4.2 we present experimental results with which we 
hope to compare our calculations. In section 4.3 we present our 
calculations and provide a discussion. Section 4.4 contains the 
conclusion . 

4.2: Experimental Results 

Table 4.1 depicts the thermal Hall energy of the DX center 
in Ga A1 As obtained by various workers. As the concentration of 

1-X X 

A1 species increases^ the band gap increases and so does At 

X == 0.45 ^ there is a cross over from the r point to the X point. 
A key result of the table is the value of AE^/Ax ^ 700 meV (Chand 
et al)^. We note in a similar work , Shubert and Ploog quote 
with a value diminished approximately by a factor of two. There 
also exists a parallel work on Sn by Basmaji using 
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magne tores istence with values of roughly the same as the ones 
quoted in Table 4.1. No clear evidence of the DX center has been 
observed for Oxygen and Carbon. The values in Table 4.1 are those 
in which no hydrostatic pressure has been applied to the sample. 
The doping is also not excessive. Further there is a small but 
perceptible chemical trend : for Si is larger than for Sn. Our 

endeavour will be to understand the trend with respect to the 
aluminium concentration in the host and also with the chemical 
characteristics of the donor species. 

Table 4.2 gives the DLTS activation energies for Group IV and 
Group VI donors. A key result of Table 4.2 is the existence of 
clear chemical trends : the Group VI donors have a const E (=s 300 
meV) while the Group IV impurities exhibit a monotonic decrease as 
we go down the column. 

4.3: Theoretical Results 

Our calculations are based on the Koster Slater formalism 
described in Chapter 2. The substitutional disorder on the cation 
site was taken into account using 
(i) The virtual crystal approximation (VGA) 

<ii) The single site Coherent Potential Approximation (CPA) as 
formulated by us. 

These are described in earlier chapters. 

We devote our attention to two systems: Ga^ .j,Al^_ 3 As <x=0.3) 

and Ga A1 As (x = 0.4). The band gap 
o . <s o . ^ 


for 


these 


systems 
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obtained by using the Vogl tight binding parameters are 1.8 eV 
and 1.95 eV respectively. These are in agreement with the 
experimental values quoted in Table 4.1. 

The defect potential on the anion site is defined as 


V 


D 


V 

ns 


- V 

os 


(4.1) 


Here is a shifted s-state ionization potential for the 

principal quantum number n of the donor species involved : n = 3,4 
and 5 for S, Se and Te respectively. These ionization potentials 
are taken from the tables of Clement! and Roetti^*. This value is 
shifted by a relatively constant value to define the zero of 
energy with respect to the vacuum and not with respect to the top 
of the valence band. V is simply the average s state anion As 

as 

potential . 


V = (i - x) E (s,a) + (x) E (s,a) (4.2) 

as GaAs ALAS 

The parameters E(s,a) are taken from Vogl et al. Calculations 
presented are accurate upto 5 meV at best. 

In figures 4.1 and 4.2 we plot the a^ anion defect level as a 
function of the inverse potential d/V^). The latter is simply the 
diagonal Green's function (at the impurity site) as the Koster 
Slater equation would demand. The VCA and the CPA plots are nearly 
identical. This result is not surprising since the relevant input 
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parameters E(Sja) of GaAs and that of AlAs differ by less than 
10%. The VGA is a good approximation under such circumstances. 

Table 4.3 depicts the calculated impurity levels for the 
Group VI elements. For Te, the value of E - E = 120 meV is in the 

C D 

same range as the experimental value of Hall energy E^ = 100 (-50) 
meV quoted in Table 4.1. The trend for increasing x is consistent 
with observation. Further the chemical trend is very weak . There 
is a reduction of only 15 meV as we go from S to Te. This is in 
conformity with the DLTS observations in Table 4.2. It is also 
borne out by Hall measurements^^. 

The more interesting observations are for the cation donor 
species. Si is by far the most extensively studied DX center. We 
may define the a^ cation defect potential in a manner analogous to 
equations (4.1) and (4.2). In figures 4.3 and 4.4 we plot the 
cation defect level as a function of inverse defect potential. The 
VGA and the GPA plots are different. This may be understood by 
examining the Vogl^*^ input parameters E(s,c) of GaAs and AlAs. 
These differ by 80% . 

Table 4.4 depicts the calculated impurity levels for the 

Group IV elements. The GPA values for Si and Sn are in good 

agreement with the experimental Hall activation energies E^ 
reported in Table 4.1. Further the average value of the slope 

obtained by us is A(E - E )/Ax ^ VOOmeV. Chand et al report an 

experimental value of AE/ix = 707meV tor SI. Moreover the 
chemical trend of a small but perceptible mcnotonic decrease in 
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with increasing Z the atomic number (as we go down Group IV ) is 
also borne out. The VGA results also exhibit this trend, but the 
absolute values of the activation energies as well as that of the 
slope - E^)/Ax are in poor agreement with the experiments. 

4.4: Conclusion 

The results of our CPA based calculations mimic the 
experimental Hall activation energies in terms of (i) Chemical 
trend as one donor species (S) substitutes another (Sn) (ii) The 
slope AE^/Ax. Both Group IV and Group VI impurities are well 
accounted for. 

Our calculations show that the DX center is a level with a^ 
symmetry. The calculations for the t^ symmetry level has also been 
done and no correlation with the experimental values have been 

13 

discerned. The work of HJalmarson and Drummond as well as 
Oshiyama and Ohnishi^* also indicate that the DX center has a^ 
symmetry. Equipped with this information about the DX state, the 
phenomenon of ppc (persistent photoconductivity) can be understood 
on the basis of symmetry arguments alone. Ga^_^Al^As is a direct 
band gap material when the aluminium concentration is less than 
45% . So in this range of A1 concentration the r point constitutes 
the minimum of the conduction band. The transition matrix element 
involves an operator which depends on position. This matrix 
element will be zero for a conduction band edge (T point) to DX 
state (a^ type). Thus the carriers will have to get excited to a 
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higher state before they drop to the DX center defect level in the 
gap, resulting in persistent photoconductivity. However when the 
A1 concentration is increased beyond 45% the X point takes over as 
the conduction band miniinuin. The X to transition is not 


forbidden and hence 

there 

should 

not 

be 

any 

persistent 

photoconductivity in 

this 

region . 

This 

is 

borne 

out by 


3 *7 

experiments ' . This does not however rule out the possibility of 

a small distortion which has little effect on E and E , and 

T e 

3 7 jL3 

results in persistent photoconductivity and large Eopi ' ' ' 


It would 

be 

interesting 

to 

explore 

the metastable 

characteristics 

of 

this defect 

by 

performing 

total energy 


calculations. The methodology using the Hellmann - Feynman theorem 
within the tight binding framework has been outlined in the second 
chapter. Pressure dependence also needs to be theoretically 
investigated. This task is reserved for future. 

Further, the DX center has been found to be EPR invisible^’. 
In view of this characteristic of the defect a negatively charged 
donor is plausible. This point is explored in the next chapter 
where we examine the band model for the DX center. 
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Experimental Hall Energy (meV) in Ga^_^Al^As as reported 
several workers. Eg is the band gap in eV. 
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TABLE : 4.2 


DLTS Activation Energies (eV) in Ga A1 As 


Author 

s 

Se 

Te 

S i 

Ge 

Sn 

Kumagai ^ 

X 

0.28 

(0.43) 

0.28 

(0.43) 

0.28 

(0.434) 

0.43 

(0.3) 

0.33 

(0.3) 

0.21 

(0.3) 

Lang-Logan 

X 

i 

0.28 

(0.4) 

0.43 

(0.4) 

- 

- 

0 . 19 
(0.4) 

Wagner^ 

- 

0.32 

- 

- 

- 

- 

Zhou^et al 

- 

- 

- 

0.44 

- 

- 

Kunzel 

L 

- 

- 

— 

0.41 
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TABLE : 4.3 

The calculated ai anion defect level in Ga A1 As. The CPA and 

1-X X 

the VGA are nearly identical, hence only the CPA results are 
shown . 
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FIGURE CAPTIONS 

Figure 4.1: Plot of inverse defect potential strength vs. the 

anion defect level in Al^, As as obtained by CPA ( ) and 

VGA ( ). 

Figure 4.2: Similar to figure 4.1 except that the A1 concentration 
in this case is 30% . 

Figure 4.3: Similar to figure 4.1 except that in this case the 
defect is at the cation site. 

Figure 4.4: Similar to figure 4.2 but with the defect at the 
cation site. 
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CHAPTER 5 : MANY ELECTRON EFFECTS IN THE 
EFFECTIVE MASS THEORY 


5.1: Introduction 

In this chapter we consider many electron effects in the 
context of Effective Mass Theory (EMT)\ We apply EMT theory for 
helium like impurities in the semiconductor matrix. This model is 
relevant to the chalcogens and some noble metal impurities in 
silicon. The discussion of helium like impurities is also relevant 
to anti-site defects in III - V semiconductors, a currently 
popular model for EL2 in GaAs^ It is further relevant to the 
study of DX center in Ga^_/l/s as this defect is presently being 
viewed by several authors as a negatively charged center . 

We present here a brief outline of the effective mass 
formalism. A perfect crystal can be described by an effective one 
electron Hamiltonian as 


Where 


H = T + V 

O ^ 

T = - h^/2m^ 


(5.1) 


is the kinetic energy and 
has the full symmetry of the 
problem is 


is the potential energy term, which 
lattice. The corresponding eigenvalue 
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Here k is the wave vector and n is the band index. 

A crystal with an isolated defect (dilute doping in 
semiconductors) can also be described in one electron 
approximation by a Hamiltonian H 


H = T + V 

Here V is the new one electron potential and it may be written as 

V = V + u 

o 

So H = H + U 

o 

Here U is the defect potential. The eigenvalue problem for an 
imperfect crystal is 


H = E 

We will be interested in those energies which are within the 
forbidden band gap. These bound state energies will be dealt 
within the framework of EMT . 

Theoretically it is possible to expand the bound state 
eigen function in terms of all Bloch functions. However, to 
simplify the task, the effective mass approximation is invoked. 
Pertaining to a given physical situation it enables us to consider 
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only a few conduction band states in terms of which the bound 
state is expressed. 

The detailed derivation of the effective mass equation in 
reciprocal space has been given by several authors^. However, the 
nature of the approximation involved in this theory can be 
understood by the following simplified model of the defect. 

Consider a substitutional donor which has one excess nuclear 
charge +e relative to the host and one excess valence electron. In 
view of the great stability of covalent bonding all electrons of 
the impurity atom except one will be engaged in bonding. This 
extra electron experiences a potential of the excess nuclear 
charge +e on the impurity. In vacuum this is exactly the situation 
of hydrogen atom. However, in a solid some modifications are 
needed for the Schrodinger equation pertaining to this problem. 
They are as follows 

(i) The nuclear charge +e polarizes the solid, leading to 
dielectric screening. Thus, the potential experienced by the 
excess electron is -e^/^r instead of -e /r . Here £ is the static 
dielectric constant of the solid. 

(ii) In the hydrogen atom the origin of energy is taken to be at 
infinite separation where the electron becomes free. Here the 
lowest energy at which the electron becomes free is the bottom of 
the conduction band E^. Hence E is chosen to be the origin of the 

energy. 

(iii) A free electron in vacuum propagates with its own mass m. 
However, in a solid 'the presence of periodic potential influences 
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this propagation. This is taken into account by introducing the 
concept of effective mass m 

Thus the Schr^dinger equation for this excess electron is 


2 m 


* 


4- e_ 
e r 


f V' = (E “E ) 
* c ' 


(5.2) 


For those impurities which have two excess electrons relative to 
the host atom the problem is similar to the helium atom. 
Consequently, the Schrodinger equation for these two excess 
electrons will be 

’ - Z'e"" I V' = (E -E ) V' (5.3) 

I 2 m ^r^ ^r^ ^ r^^ J 

Here Z' is the excess nuclear charge. Note that for negatively 
charged center (as is the case of DX ) the excess nuclear charge 
will not be equal to the excess number of electrons. The above 
equations are valid for one band with one extremum. 

The effective mass equation for a band with several extrema 
was first derived by Twose (reported by Fritzsche) . If one 
neglects the intervalley mixing, then each minima can be 
considered as being independent of others. One obtains the usual 
effective mass equation for each minima separately with the 
effective mass tensor corresponding to that particular minima. If 
however, all the minima are equivalent then the effective mass 
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tensors corresponding to all the extrema will be identical leading 
to a degenerate level in the gap. The degeneracy of the level is 
the same as the number of equivalent minima. This degeneracy is 
lifted by incorporating the intervalley mixing. There are several 
ways of incorporating this effect. One such approach is the 
conventional first order perturbation theory described by Bourgoin 

<S V 

and Lannoo . Resca and Resta have shown that in the case of 
intervalley mixing the effective Schrodinger equation of the 
envelope wave function associated with the donor is the potential 
V(r) modulated by a spatial renormalization factor j<^(0,r)l^. We 
will be neglecting the intervalley mixing effect in all our 
calculations . 

In section 5.2 we present our formalism for treating helium 
like impurities. Section 5.3 contains results. Section 5.4 
constitutes the conclusion. 

5.2: Fornalisn 

A great deal of work has been done on hydrogenic EMT . 
However, a corresponding effective mass theory for helium like 
impurities as defined in section 5.1, remains to be explored. 

B 

The helium like impurities were first treated by Glodeanu in 
the framework of EMT. He employed a fitting procedure for Z, the 
excess nuclear charge to account for the chemical shift observed 
experimentally. The defect level energy was obtained by a 
variational treatment of the Hamiltonian. Our endeavour will be to 
account for the chemical shift in this many body problem. We will 
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consider the impurity to be surrounded by a sphere of cavity 
radius r^. The influence of the impurity within the cavity is 
unscreened by the medium, hence inside the cavity the electron 
experiences the unscreened Coulomb potential of the excess nuclear 
charge. Beyond this cavity the medium dominates and the potential 
due to excess nuclear charge is only partially experienced by the 
electron . 

The potential energy function can be divided into the 
following three regions: 


Vj = - Ze 


Ze + 2Ze' 


£ r 
o 1 


£ r 

O 2 


£ r 
o o 


1 - £o/£ 


(5.4) 


for 


r & r^ < r. 
12 0 


II 


= zf -Ze^ - Z^ + Z^ [ 1 - -eo/'c]] 


(5.5) 


£ X ^ 'X^ 

O 1 2 O O 


for 


r < r St r, > r _ 

1 O 2 O 


III 


- Ze" - Ze 

£ r. 


£ X. 


(5.6) 


^ ^ "2 ^ 


for 
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The above equations give the potential experienced by the 
electron excluding the interaction term, which will be dealt with 
later. Here and c are the dielectric constant for the vacuum 
and the medium respectively. The third term in and has been 
added specifically to make the potential continuous at the 
boundaries of the three regions. Note that this term is a 
constant. The interaction term (where ^ or 

12 O 

depending on the values of r^ and r^ ) can not be made continuous 
by a mere addition of a constant term. However, for a given 
probability charge distribution Ivl^e of a particular electron the 
potential experienced by the other electron can be made 
continuous. 

The above problem was solved using two types of trial 
functions. The results in both the cases were found to be in fair 
agreement with each other. 

In the first case the trial function was chosen to be the 
product of two hydrogen atom wave functions 

V (r^,r 2 ) = l/n (n/a)^ exp(-Tjrya) exp(-77r2/a) (5.7) 

Here t? is the variational parameter. The electron electron 
interaction was chosen to be e when both r^ and r^ were 

less than ro and eV'Cr^^ otherwise. The potential experienced by 
one electron due to the probability distribution of the form 
(Va)®'^'^exp(- 77 r/a) of the other electron can be made 
continuous at the boundary of the cavity by the addition of a term 
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dependent on 'r?' . 

The average value of the Hamiltonian in the state ) is 

evaluated analytically. Newton Raphson method is employed next, to 
find the value of the parameter t) at which the energy is minimum. 

Figure 5.1 exhibits a typical plot of ground state binding 
energy per electron as a function of the cavity radius r . Here ^ 

m 

and m are taken to be 10 and 0.5 respectively. The energy level 
is shallow for small values of r^. At r^ equal to 1.4 a.u there is 
a crossover from shallow to deep level. This shallow to deep 
instability is often encountered in experimental situations (see 
section 5.3), It may be understood simply: when r^ the cavity 

radius is relatively large, the electron is in a free atom 
environment. When r^ is small, the semiconductor matrix dominates. 
The crossover point r^^ is found to be a monotonic increasing 
function of £'and monotonic decreasing function of m . The 
effective mass of GaAs is less than that of AlAs, and the static 

dielectric constant of GaAs i.s more than that of AlAs. This 

results in lowering of ^oc Soss towards higher 

concentration of A1 in Ga^_^Al^As. 

The behaviour of energy with respect to the cavity radius r^ 
is similar in spirit to that observed by Bourgoin and Mauger in 
the context of hydrogenic EMT with a cavity. The above choice of 
trial function has been employed to arrive at the results of 

section 5.3. 

In the second case the trial function is chosen to be the 

standard variational solutions for 


linear combination of two 
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helium atom immersed in vacuum and in the medium. They are as 
follows . 



= l/rr (h/a)^ exp(-r)^r^/a) exp( -rj^r^/a) 

(5.8) 


= 1/*^ i'n /3.) exp(-r),r /a) exp(-7? r /a) 

(5.9) 


where = (Z-5/16) 



and ri = (Z-5/16)m*£ /m £ 

^ O 



The trial function is 

^'(r^^r^) = (5.10) 

This choice of trial function is similar in spirit to Reiss^°. He 
has chosen a similar trial function to solve the problem of 
hydrogenic EMT with a cavity. 

In this case the expansion coefficient is the variational 

parameter. When r^ is zero the situation is similar to helium atom 

immersed in the medium with no cavity. In this case is found to 

be zero and C is unity. Similarly when r is infinity implying 
2 

that there is no effect of the medium, is found to be zero and 
is unity. For an intermediate value of r^ both terms have non 
zero coefficients. 

In the second method the potential experienced by one 
e-lectron due to the electron electron interaction term cannot be 
made continuous for both 4^^ and <P^ simultaneously by adding the 
same constant. Hence, the dielectric constant was taken to be ^ 
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throughout for this particular derivation. The expectation value 

of the Hamiltonian in the two states 4 > and 0 can be obtained 

1 2 

once again analytically. 

Def ining 

«« = < ■#'J H 

= < 0 J H 

H,, = < 0.1 H 10.> 

= < 0J H |0.> 


The ground state energy is obtained from the solution of the 
following equation. 


H - E 
11 


H - A E 
21 2 1 


(5,11) 


Hi2 - ^lE H22 - E 


Figure 5.2 shows a typical plot for the ground state energy per 
electron vs. the cavity radius r^ for this particular derivation. 
The result agrees with that of case 1, implying thereby that an 
approximate treatment of the interaction term is fairly justified 
in this context. Small changes in this term do not alter the 


result drastically. 
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5.3: Applications 

5.3a: Chalcogen Donors in Silicon 

The chalcogen donor level in silicon is six fold degenerate. 
Since the chalcogen occupies a substitutional site, the overall 
symmetry is tetrahedral (T^). The symmetry is broken by the 
central cell potential and is also known as valley-orbit 
splitting. The six fold degeneracy is thus lifted to a singlet a^, 
a triplet t^ and a doublet e. The singlet a^ is the lowest. Spin 
orbit effects further split the t level to r (doubly) degenerate 
and (non - degenerate). Stray fields are responsible for 

lifting the degeneracy of the doublet states e and r . 

8 

The experimental situation has been exhaustively reviewed by 
Grimmeiss and Janzen^^. The chalcogen species S, Se, Te have been 
known to introduce two distinct levels. The deeper level known as 
the A-center is now identified with the ionized impurity atom, 
e.g. S"^, Se"^ or Te^. The shallower ones, labelled B-center 
correspond to the neutral species S , Se , Te . 

Our calculations will focus on the neutral B-center. Further, 
only the ground state energy is considered. Table 5.1 depicts the 
current experimental status of the chalcogen species in Si 
(excluding Oxygen). For completeness, we also display the energies 
attributed to chalcogen related defects. Our interest is confined 
to the second column of Table 5.1. The calculation of excited 
state energies, for example (Is - ai) ^ (Is tz) within 

effective mass theory can be performed."^ We would, however, first 
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like to examine the ground state energies within our proposed 
cavity model described in section 5.2. 

In figure 5.3 we plot the double donor level of Si as a 
function of cavity radius. The cavity radius has been expressed in 
units of the nearest neighbour distance r for silicon (r “ 4.43 
a.u). Figure 5.3 shows that there is very little change in the 
position of the defect level as the cavity radius is increased 
upto approximately 0.665 times r . This is the region where the 

nn 

medium dominates. Beyond this as one increases the cavity radius, 
the position of the defect level shifts very rapidly, the energy 
level curve shows a very steep change and the level becomes 
deeper . 

Our endeavour will be to explain the results of column 2 of 
table 5.1 within the framework of EMT . We need an appropriate 
value of the cavity radius for different chemical species that are 
depicted in table 5.2. One method of determining the cavity radii 
is to employ the innermost orbitals of the doubly ionized atomic 
species. We define an average size of the doubly ionized species 
as 


<r > 

av 


Z < r . > n . / ^n . 

t t t t b 


Where, for 

sulphur is equal to Is, 2s, 2p, 

3s , 

3p ^ 

for 

selenium (Se^^) 'i' is Is, 2s, 2p, 3s, 3p, 3d, 4s, 

4p; 

and 

for 

Tellurium 

(Te^^) 'i' is Is, 2s, 2p, 3s, 3p, 3d, 4s, 

4p, 

4d, 

5s , 

5p . These 

radii have n, as 1 for s, 3 for p and 5 for 

d. 

For 

the 

outermost 

valence p-orbital is 1. These radii are 

tabulated 

in 
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table 5.2. Hote that these are distinct from the neutral atomic 
radii. The latter increases with atomic number Z and primarily 
determines valence based solid state properties. The core radii in 
tabl© 5.2 are typical of the atom and not of the solid. 

The agreement between the calculated and the observed values 
is fair. For impurities like sulphur, selenium and tellurium, 
whose core radii are found to be larger than the crossover cavity 
radius the error range for the defect energy is found to be 
large. Due to the steep rise of the defect energy as a function of 
r^, a difference of. 0.001 a.u in the value of <r > can give rise 
to a difference of about 0.05 eV in the defect energy. Thus a 
rounding error for <r > even at the third decimal place gives 

av 

4. 

rise to an error of about - 0.05 eV in energy. In this context the 
results reported for Magnesium are fairly accurate. The energy 
curve at the core radius of magnesium is fairly steady. Thus minor 
differences in determining the core radius do not get reflected in 
the defect level energy. 

5.3b: The DX Center: A Negatively Charged Donor 

As mentioned in the concluding section of chapter 4 the DX 
center defect which is found in Ga^_^Al^As is EPR invisible. In 
view of this characteristic it is conjectured that this defect is 
due to a negatively charged donor^ A negatively charged donor 
will have two electrons in excess to the host atom which it 
substitutes. Hence, the analysis of section 5.2 is applicable for 
an investigation of this situation. 
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We present preliminary results. Figure 5.4 depicts the donor 
level as the function of cavity radius for Ga A1 As The r 

0.7 0.3 

band related donor level is uniformly shallow. The X band donor 
level is the deepest, being deeper than the L band. There appears 
to be no unique way of fixing the cavity radius for the negatively 
charged center. If for the sake of argument we fix the cavity 
radius at a value of approximately half the nearest neighbour 
distance (say r^ % 3a. u) then figure 5.4 predicts a DX level tied 
to the L band and positioned approximately at 300 meV below the 
conduction band. There is a shallow level tied to the r band. The 
defect level tied to the X band is inside the valence band. The 
same value of r^ for GaAs (in figure 5.5) leads to two shallow 
levels which track r and L band respectively. The level due to X 
band is nevertheless inside the valence band. 

In this context it is worthwhile quoting a similar work 
based on EMT by done Burgoin and Mauger'* for single electron 
donors. They have incorporated the multivalley mixing effect of 
several equivalent minima for X as well as L bands. Due to this 
intervalley mixing effect the L derived state with a^ symmetry is 
found to be deeper than the X derived state of the same symmetry. 

We have shown above that a simple cavity model without 
incorporating the correction due to intervalley mixing can explain 
the L band tracking of the DX state. Further, the conjecture of it 
being a negative center, put forward to explain 
invisibility is also borne out by our calculations. 


the EPR 
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5.3c: The EL2 Center: an Anti-site Defect 

EL2 is the commonest native deep defect in GaAs . It is 
called the main electron trap. It is conjectured that it is an 
anti-site defect namely, arsenic on the cation gallium site 
(ASj^g^). Calculations using the Koster Slater equation have 
reported this defect to be deep. 

It is clear, however, that ASq^^ is a substitutional helium 
like donor. Hence it is worthwhile applying the analysis of 
section 5.2 to this problem. Figure 5.4 represents the plot of a 
double donor level with respect to the cavity radius r^ . We note 
that the nearest neighbour distance in GaAs is r ^4; 5.6 a.u. A 

nn 

cavity radius of about r^ % 4r yields a deep level at the mid 
gap. Thus our calculations show that an effective mass approach to 
the problem can possibly explain the deep level characteristic of 
the EL2 level. 

5.3d: Double Donor in Ga^ A1 As 

In this context it is interesting to observe that as we 
approach AlAs the double donor level becomes deep more readily 
(e.g at a lesser value of ro). This can be seen by examining the 
plot of Ga A1 As in figure 5.7 and comparing it with the 
plot of GaAs in figure 5.6. In figure 5.7 the crossover from 
shallow to deep defect level starts when the cavity radius is 
around 7 a.u whereas for pure GaAs it occurs at about 11.4 a.u. 

Fixing the cavity radius at a particular value one observes 
that the donor level gets deeper from the conduction band as the 
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aluminium concentration is increased from 0.0 to 0.45 (see figure 
5.8). For aluminium concentration higher than 0.45 the X band 
takes over as the minimum of the conduction band. Thus beyond x = 
0.45 an analysis relevant to the X band needs to be done to view 
the trend of the double donor defect level with aluminium 
concentration. For aluminium concentration less than 45% we see 
that the same impurity (characterized by a particular cavity 
radius) will give rise to an increasingly deep level in the 
environment of increased aluminium concentration. We believe that 
similar trends will' be observed for hydrogenic donor impurities. 

5.4: Conclusion 

In this chapter we have proposed a simple cavity model 
Effective Mass Theory for helium atom and have applied it to 
several problems of interest. 

We have studied the chalcogen impurities in silicon and have 
obtained results which agree well with experimental values. 
Further the DX center, which is currently being viewed as a 
negatively charged center due to its EPR invisibility, has also 
been studied using this methodology. Our results do not rule out 
such a possibility. However, a firm statement in favour of such a 
conjecture cannot be made since the cavity radius in this case 

remains an empirical parameter. 

The anti-site defect in GaAs has also been treated. We have 
found that for a reasonable value of the cavity radius, the 
theoretical results agree with the experimental values for the EL2 
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defect center, suggesting thereby that it could be an anti-site 
defect. We have also shown that the double donor level in 
Ga A1 As will become increasingly deep with increased A1 

jL“" X X 

concentration. Similar trends can also be discerned for single 
donor impurities in Ga^_^Al^As. 

Thus, even a simple theoretical model of helium atom within a 
cavity can be successfully employed for several semiconductor 


problems . 
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TABLE 5.1 


Ground state binding energies in eV of chalcogen isolated 
substitutional atoms (D), chalcogen pairs and chalcogen 

related defects 


DEFECT 

o'* 

d"" 


"VI 

c ‘ 







O . 11 


S 

0.32 

0.61 

0.19 

0.37 

O . 16 


1 

Se 

1 

0.31 

0.59 

0.21 

0.39 

0.12 

0.21 

Te 

0.20 

0.41 
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TABLE 5.2 


The core radii of the chalcogen impurities and of Mg in 

atomic units. The calculated binding energies E are read off from 

b 

figure 5.3 associating <r> with cavity radius r The 

<3'V 0 ‘ 

experimental binding energies are taken from Table 5.1 and 
reference cited therein. For Mg, the experimental value is the one 
reported by Milnes. The binding energies are in eV. 


IMPURITY 

CORE 

RADIUS 

<r> 

av 

THEORETICAL 

EXPERIMENTAL 

S 

3.037 

0.33 

(± 0.04) 

0.32 

Se 

3.016 

0.29 

(± 0.04) 

0.31 

Te 

2.980 

0.188 

0.20 

Mg 

2.527 

0.9092 

0.11 
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FIGURE CAPTIONS 

Figure 5.1: A typical plot of ground state binding energy per 
electron vs. the cavity radius r^ . The dielectric constant and the 
effective mass of the electron is taken as 10^: and 0 . Sin 

O * e 

respectively. This particular plot is for the product trial 
function described by equation 5.7. Note the steep rise of the 
binding energy for r^ > 1.4 a.u. 

Figure 5.2: Similar to figure 5.1 except that the trial function 

in this case is described by a linear combination as described in 

equation 5.10. Note the steep rise in the binding energy for 

r 5: 1 . 4 a . u . 
o 

Figure 5.3: Plot of double donor level in silicon as a function of 
cavity radius. The cavity radius has been expressed in units of 
nearest neighbour distance of silicon ( r^^ = 4.43 a.u ). Note the 
sharp rise in the binding energy at around 0.665 r^^- 

Figure 5.4: Plot of a negatively charged donor level in Ga^ ,j,Al^_ ^As 
as a function of cavity radius r^ . The crossover behaviour occurs 
at different values of r^ for levels tied to the X and L bands. 
There is no crossover behaviour for the r band for reasonable 

values of r^. 

o 


Figure 5.5: 


Similar to figure 5.4. The 


system in this case is 
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GaAs . The r band result remains parallel and is very close to the 
X axis. 

Figure 5.6: Plot of binding energy per electron for a double donor 
defect level in GaAs. The level is for the r band. The crossover 
occurs for a large value of the cavity radius. 

Figure 5.7: Similar to figure 5.6 for ^^^Al^ ^^^As- 

Figure 5.8: The shaded lines in the figure represent the 
conduction and the valence bands as the A1 concentration is 
increased from 0% to 45%. The continuous line shows the plot of a 
typical defect level characterized by the cavity radius r^ 

10 a.u. Note that the level deepens as the A1 concentration is 

increased . 
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FIGURE 5.8 
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CHAPTER 6: EFFECT OF DISORDER IN 
THE ANALYSIS OF EXPERIMENTS 

6.1: Introduction 

Deep level transient spectroscopy (DLTS)'*’ has been widely 
used over the past ten years as a tool to study deep levels in 
semiconductors. The context of the original proposal was that of 
defect states with well defined energy levels in crystalline host 
materials. In the past few years there has been considerable work 
on systems where either the host or the defect is disordered. In 
such cases the DLTS spectrum may be broadened and the underlying 
capacitance or current transient may be non-exponential. 
Nevertheless, in most cases, the standard DLTS expression derived 
for exponential transients is used to analyse the data. In this 
chapter we analyse the effect of gaussian broadening of a deep 
level energy spectrum on the DLTS signal. We show that in all 
cases except extreme broadening the standard DLTS analysis 
introduces negligible error. We shall also discuss the effect of 
disorder on techniques complementary to DLTS, in particular Paired 
Temperature Spectroscopy (PATS). Thr interplay between 
metastability and disorder in the interpretation of the above 
mentioned experimental techniques is also discussed. 

A semiconductor alloy, e.g. Ga^_/l/s, may be considered to 
be a substitutionally disordered system. Ga^_/l/s has a basic 
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zinc-blende structure where the cation site is occupied by Ga or 

A1 atoms. The occupancy is essentially random if the probability 

that a given cation site is occupied by a Ga(Al) atom is l-"x(x). 

However, fluctuations from the mean occupancy are possible. These 

fluctuations arise due to a variety of causes. Some of them are: 

(i) the Pauling electronegativity difference between Ga and As is 

0.4 while the difference between A1 and As is 0.5, thus favouring 

A1 as a neighbour of As. (ii) local inhomogeneities such as the 

presence of voids and anti-site defects. These can result in local 

/ 

deviations from the equilibrium concentration. A defect level in 
the fundamental gap of the semiconductor may be identified with 
its corresponding defect potential via the Koster-Slater equation 
as mentioned in earlier chapters. For the reasons discussed above, 
the defect potential in a semiconductor alloy has fluctuations 
which cause the delta function like gap level to be broadened. 

There exists a growing body of experimental data on gap 
levels in semiconductor alloys. Some of these indicate the 

existence of potential fluctuations and consequent broadening of 
the gap level. In the present chapter we shall discuss this 

broadened gap level in the context of DLTS* and other 

complementary techniques. Let us mention briefly experimental work 
other than DLTS, which indicates the existence of broadening. 
Mooney and co-workers^ have studied the kinetics of charge 

trapping in semiconductor alloys and have hypothesised a 
distribution of energy levels and capture cross sections. Raman 
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spectroscopy, too, suggests the importance of alloy potential 
fluctuations . Samuelson et al^ have attributed the splitting of 
photoluminescence lines of Cu in GaAs P to alloy potential 

l-X X 

fluctuations . 

Over the past few years there has been a growing body of 
literature on DLTS studies of disordered systems. Omling et al^ 
have studied the broadened DLTS spectra in GaAs P . Murawala**, 

7 

Singh and co-workers have discussed the broadening of the E3- 
like level in GaAs^_^Sb^. Kanieswka and Kanieswka* have similarly 
studied the main electron trap in GaAs P . Lang^ and co-workers 

1-X X 

have studied the gap states in amorphous hydrogenated silicon 
(a-Si:H) where the disorder-induced broadening is expected to be 
more pronounced than in semiconductor alloys. 

Linder ideal conditions, the capacitance transient C(t) of 
deep level is a simple exponential^'"^*^ 

C(t) = Co exp(-t/T) (6.1a) 

1/T = AT^ exp(-E/kT) (6.1b) 

Here E is the activation energy and A is related to the capture 
cross section of the deep level under consideration. T is the 
temperature, k the Boltzmann constant and r is referred to as the 
time constant. For the case of a broadened defect level the most 

7 

general expression has been considered by Singh 
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C(t) = 


Co 


2naa± 


dE 

o 


00 

dA X 
o 


exp 


(E -Eo) (A - Ao) 


2 a 


2 


T(A,E) 


In the above expression the effect of the potential fluctuations 
on both the activation energy and the capture cross section has 
been taken into account. The broadening assumed is gaussian with a 
mean at E and A , and the half width at half maximum (HWHM) of 

O O 

approximately cr and respectively. A detailed examination^ has 
indicated that the broadening in the quantity related to the 
capture cross section may be ignored. We then obtain 


dE exp{- ) (6.2) 

5-B 

This expression has been widely employed in recent years 
In this chapter we first demonstrate how expression (6.2) implies 
a long lived and non-exponential capacitance transient. We 
demonstrate this by independent analytical and numerical work. 
Next, and more importantly, we show under what conditions the 
standard interpretation of DLTS data is valid for a broadened gap 
level. At this juncture we will define the term 'standard 
interpretation'. The time constant for the apparent peak of a DLTS 
signal was shown by Lang to be^ 


C(t) = 


Co 


r r>i 


2 , 1/2 


ln(t2/ti ) 
tz - tl 


(6.3) 



6-5 


This expression contains no information on the broadening or 
smearing factor cr , The use of this expression for extracting the 
defect parameters (activation energy and the pre-exponent A) will 
be deemed the standard interpretation. 

We identify the natural scales in the problem; 


£ = Eo/kT 
s = <y/Eo 
b = a/\Cl 


(6.4a) 

(6.4b) 

(6.4c) 


Unless the level is very shallow, we have £ >> 1 and this scale 
shall not concern us henceforth. If the potential fluctuations are 
small, then s and b in equation (6.4) will be small. We shall call 
the case { b<l, s<0.1 } the weak disorder case. The regime {b2:l, 
s^O . 1 } will be deemed the strong disorder case. We predict 
significantly different behaviour for these two cases. 

A complementary technique for defect analysis based on 
capacitance transients was first proposed by Singh et al*^. It is 
Paired Temperature Spectroscopy (PATS). In PATS the signal is 
defined as a function of time and its value is the difference 
between the transients at two temperatures. The major advantage of 
this technique is that it does not need a continuous thermal scan; 
two to three fixed temperature baths are sufficient for this 
purpose. The fixed temperature permits the use of signal averaging 
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and an enhanced 
the PATS signal 


signal to noise ratio can be achieved, 
occurs at a time 


The peak of 


ln(T /t ) 

1 2 .' 

OTFTr-XTTT) 


In section 6.2 we treat the case of weak disorder {b<l, s< 
0.01}. We show that the capacitance transient is weakly 
non-exponential. We validate the standard DLTS analysis based on 
equat ion ( 6 . 3 ) . In section 6.3, where we treat the strong disorder 
case {b ^ 1, s ^ 0.1}, the transient capacitance is highly 
non-exponential, akin to glassy and spin-glass relaxation. 
Further, we demonstrate that the standard DLTS analysis will 
systematically overestimate the activation energy in cases of 
severe broadening. In section 6.4 we apply our experience of the 
previous two sections to two realistic cases: a-Si:H previously 
studied by Lang et al*^ and an E3 like level in GaAs Sb 

O . B<S O . 14 

previously studied by Murawala et al*^. 

As noted in chapter 4 the DX center as well as several other 
defects are metastable. In section 6.5 we present a novel analysis 
of the DLTS peak heights for a metastable species in the presence 
of disorder. Section 6.6 encapsulates our conclusions on related 


phenomena . 
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6.2: Analysis for Weak Disorder 

We first consider the weak disorder case {b = <y/kT < 1, cy/Eo 
< 0.1} in DLTS. As cr tends to zero the distribution approaches a 
delta function 6(E - Eo). We are then justified in expanding the 
exponential term in the integrand of equation (6.2) 


exp(-E/kT) = exp(-Eo/kT) exp(-(E-Eo)/kT) 

^ exp(-Eo/kT)x[i-(E-Eo)/kT + (E-Eo)V2(kT)’'] 


Def in ing 

1/T^ = AT^ exp(-Eo/kT) (6.5) 

C(t) = Co exp(-t/To) (6.6) 


equation (6.2) becomes 


Eg 


C(t) % C.(t) 


dE 


/ o 2,1X2 

o (2n<y ) 


exp 




(E-Eo) 


t(E-Eo)o- 
kT To 


■}] 


We now apply the stationary point formula to the above 
integral. Given 


1(a) = J exp( -a<;i;(E ) )dE 
the formula^^ states that for large ot 
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1(a) % (2n)^^^exp( - a0(Ep))/| - a0''(Ep) 

where Ep represents the minima of <^(E). Application of this 
prescription yields 


C(t) 


Co exp I - + 0-5(cyt/kTTo) 

L 1 +<o'/kT)^t/To 



+ (0'/kT)^t/To] 


(6.7) 


Equation (6.7) is our modified formula for the case of weak 
disorder. Note that for the case a tends to zero we regain the 
exponential expression (6.6) for the ideal transient. 

For a DLTS analysis we construct the signal 


S(T) = (C(ti) - C(t2))/Co 


( 6 . 8 ) 


and differentiate it with respect to the temperature to obtain 
(see Appendix) the equivalent of equation (6.3) 

1 _ ln( tz/ tl ) ^0 

(t2-tl)[l +(a/kT )^] 

m 

The subscript m denotes the maximum DLTS signal. Expression (6.9) 
reduces smoothly to the standard DLTS expression (6.3) as o'/kT^ 
tends to zero. The implication is that the standard DLTS analysis 
may be in slight error and yield an incorrect activation energy. 
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For disordered systems the analysis of PATS is similar to 
that of DLTS. The maxima of PATS signal for weak disorder occurs 
at a time 

= ln(T^/Tp/[l/T^(l + (a/kT^)^) - 1/T^(1 + (cr/kT^^)] 

We see that the expression is identical to that of DLTS except for 
the interchange between the time and the inverse of the time 
constant. The inferences drawn in this case are the same as that 
of DLTS. 

Figure 6.1 depicts the capacitance transient as calculated 
using the exact expression(6 . 2 ) , our approximate equation(6 . 7 ) and 
the pure exponential equation(6. 1) . The three curves are nearly 
identical, though our approximate expression (6.7) and the exact 

(6.2) are in closer agreement. Further, the exponential expression 
decays faster than the other two, indicating the slowly decaying 
(non-exponential) character of the broadened gap level. 

Figure 6.2 depicts the DLTS plot for the gap level discussed 
in section 6.1. The maxima of the signals are at identical 
temperatures for the three cases discussed in the previous 
paragraph. This implies that the modified expression (6.9) has 
little influence on the determination of the activation energy. 
For weak disorder the standard DLTS analysis based on equation 

(6.3) would yield accurate values of the activation energy and the 
pre-exponential factor. Thus we conclude that the 


standard 
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analysis is robust with respect to small disorder. We have indeed 
this for several values of the parameters b and o'/Eo. 

6.3: Analysis for Strong Disorder 

We now consider the case of strong disorder {b = o'/kT > 1, s 
- (y/Eo > 0.1}. Defining x = (E-Eo)/©', and using the definitions 
(6.5) and (6.4), the capacitance transient of equation (6.2) may 
be cast in the form 

00 

C(t) •= Co J exp(-f(x))dx/(2n)^^^ (6.10) 

-Eo/O- 

where 

f(x) = xV2 + (t exp(-bx))/To (6.11) 

On account of the nature of the gaussian, the upper limit E^ in 
the equation (6.10) has been replaced by infinity. This entails 
only negligible error. We next apply the stationary point formula 
in a manner similar to the previous section. Differentiating, we 
get 

f'(x) = X - bt exp( -bx)/ro 

For large b the solution of f'(xo) = 0 is given by the recursive 


expression 
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Xo 


In 


tb /To 


In 


tb^/To 


ln(tb /To 


) 


The above expression holds under the condition tb^/To > e. We 
shall use a truncated form of the above expression, i.e. stop the 
recursion after the first denominator. A series of algebraic 
manipulations then yields 


C(t) % Co exp 




tb^/To 


In tb /To 




+ ln(tb /To) 




1^1 + ln(tbV'J'<=) 


( 6 . 12 ) 


Equation (6.12) is a radical departure from the Debye 
relaxation (pure exponential) form. If the first term in the 
exponent is ignored, it reduces to 

C(t) « l/t“ (6.13) 

where a = 1/b^, a form of non-exponential relaxation currently 
popular in spin-glass and related areas. 

Figure 6.3 depicts a typical capacitance transient 
encountered in a DLTS analysis of a-Si:H^. Here b = o'/kT is large 
(% 6). The transient is long lived and distinctly non-exponential. 
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The approximate expression (6.12) agrees moderately well with the 
exact one. 

We have been unable to obtain a closed-form expression for 
the DLTS time constant extrema analogous to equation (6.3) or 
(6.9). The difficulty is due to the complicated character of 
equation (6.12). Hence, for further analysis of the strong 
disorder problem, we have to rely on the exact expression (6.2). 
The DLTS extrema can be located by differentiating equation (6.2) 


as 

ar 


= J dE 'exp [l-.(E-Eo)V2o'^] x [exp( -ti/r )/T^-exp( -tz/r )/t^] 


= 0 


(6.14) 


and employing a non-linear numerical procedure such as 
Newton-Raphson to extract the three parameters E^, A and o-. 

A question of central importance is the error made if the 
standard analysis based on equation (6.3) is applied to the DLTS 
data for the strong disorder case. The analysis of the previous 
section indicates that the error made is almost neglible for the 
weak disorder case. We consider a typical deep level centred at 
= 0.5eV, and A = 10%■V^ and characterised by an increasing 

broadening c. -- 0 . 1 , 0.25 and 0.4eV, respectively. Using a rate 

window tz/ti = 4, DLTS signals employing the exact broadened 
expression (6.2) were generated. From the maxima of these signals, 

w using the standard expression 

T , we construct Arrhenius plots usi g 

m 
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(6.3). Figures 6.4 - 6.6 depict Arrhenius plots for the three 
distinct broadening considered {<y - 0.1, 0.25 and 0.4eV 

respectively). An examination of figures 6.4 - 6.6 leads us to 

conclude (i) there are no deviations from linearity over the four 
time decades considered in the DLTS simulations (O.lms to Is). 
Hence an examination of the Arrhenius plots provides little clue 
to the broadening of the gap level, (ii) The extracted gap level 
E ' is invariably greater than E , the actual peak in broadened 

O O 

density of states. In fact as figure 6.6 depicts, for o* = 0.4eV, 

Eo' = 0.72eV, while Eo = 0.5eV, a 50% deviation! (iii)The 

extracted value of the pre-exponent A' is also larger than A. 
However no systematic trends were found for A. 

Figure 6.7 displays the extracted value Eo' as a function of 
<//Eo. It reinforces and encapsulates the conclusion based on 
figures 6.4 - 6.6. For o’/Eo > 0.1 there is a systematic 

overestimation of the activation energy. For a/Y.o <0.1 this error 
is negligible. 

An analysis for strong disorder in PATS has been carried by 

19 

Singh et al . 

6.4; Simulation on Real Systems 

We shall apply our analysis to two experimental situations: 

(i) a-Si;H (ii) an E3 like level in 6-aAs^ ^^Sb^ 

For the a-Si:H case^ the disorder is large and a priori we 
be large. Indeed, as figure 6.8 shows, the fit to the 


expect O' to 
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experimental DLTS signal is accomplished by using o' = 0.175eV, 

Eo = 0.82eV and A = 4x10* s ^ in equation (6.2). This implies 
that o'/kT^ ~ 6 and o'/'Eo ^0.2. This fit was possible only after 
some trial and error. The corresponding Arrhenius plot is depicted 
in figure 6.9. Note that it is linear over five time decades. The 
value of the activation energy extracted from the slope of the 
data presented in figure 6.9 is Eo' = 0.87 eV . Recall that the fit 
in figure 6.8 used the input value of Eo = 0.82eV. Thus the 

standard DLTS analysis overestimates the activation energy in this 
case by O.OSeV. Lang et al^ have done numerical analysis for the 
case of a-Si;H, taking into account both large trap densities and 
broadened energy spectra, and find similar results. Thus, our main 
thesis encapsulated in figure 6.7 that the standard DLTS 
overestimates the activation energy in the limit of severe 

broadening, is borne out. 

The next case we consider is the E3 like level in 

GaAs Sb ■ This has been subjected to a DLTS study by 

O , 8<5 O . 14 

Murauala and ac.-»orkers°. A close fit to the experimental DLTS 
plot is obtained by simulating equation (6.2) uith Eo = 0.23eV, A 

= 2. SlxlO-s^'r'. and <^ = O.OSSeV. This corresponds to o'/kT™ = 4 

and o/Eo =0.24. »e have not reproduced the DLTS signal fit here 

since it is analogous to figure 6.8, Hooever, »e do sho» the 

Arrhenius plot in figure 6.10. The plot is again linear over five 

time decades . The extracted activation energy is 0.25 'eV. in 

slight excess of the peak in the broadened gap level (Eo = 
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0.23eV) . 


6.5: Disorder and Metastability 

Levinson has reported an analysis of the DLTS peak heights 
for a defect oscillating between two metastable states . In this 
section we examine the effect of disorder and present some new 
resu Its . 

In figure 6.11 we present an idealized model for a 
metastable defect set forth by Levinson^"^- The defect transforms 
from the state B to the state A during the forward (or zero ) bias 
pulse of width t . It transforms from A to B during the reverse 
bias time of width t^ . Let the fraction of the species in states 
A and B at time t be f^(t) and fg(t) respectively. The 
normalization 

f^(t) + fg(t) = 1 (6.15) 

holds at each instance of time . Consider the time period 0 < t < t^ 
The transformation A — > B is described by a simple exponential in 
the absence of disorder 




(6.16) 


where t . is the relevant time constant, 
ab 

effect of disorder. Disorder dictates 


Now let us examine the 
that the sharp delta 
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function like level characteri 2 ing a defect is broadened . This is 
indicated by the shaded regions around the A and B states in 
figure 6.11 . We assume without loss of generality that the 
relevant density of states (DOS) is rectangular. 

n^(E) = 6(E - E^) 

disorder. , ^ ^ v 

A 0<E<E^ (6.17a) 

and similarly for ng(E) . Normalization of DOS implies . 

A E^ = 1 (6.17b) 

We wish to state that the results of this section are independent 
of this particular choice of a rectangular DOS. Equation (6.16) 
must be averaged over disorder . Hence 

fA(ti) = J n(E) f^(t.) dE 

= f° J n(E) dE (6.18) 

The time constant is given by the barrier height formula 

l/%b = ^ab 

We define a new variable x 
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tj^/T(E) = exp(-(E-x)/kT) (6.19) 

This facilitates the reduction of the exponential in equation 
(6.18) to a step function 

exp(-t^/T^^) ~ e(E - X) 

where 0(E - x) = 1 for E > x and zero otherwise . Thus equation 
(6.18) reduces to 


f^(ti) = f° Jn(E) 0(E -x)dE 
= f° J n(E) dE 

X 

3 f^ A (E^-x) (6.20) 


where the last step employs equation (6.17). Using equation (6.19) 
for X we obtain 


e 


AB 


A (E^ - kT in 
1 - kT A In C^ab^i^ 


<6. 21a) 

(6.21b) 


Consider the time period t^^ < t < t^+t^ when B ►A . An 


analysis 


similar to above yields. 
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l^(t. + tp) = 1 - (6.22a) 

= 1 - kT B In (6.22b) 

The equations (6.21a) and (6.22a) possess a form similar to 
Levinson except that and ©g^ are each different in content due 

to the effect of disorder . 

Proceeding over further cycles of reverse and forward biases 
we obtain a generalized form 

_ '00 

f^ = Cl + C2kT In(v^^tp) (6.23) 

In equation (6.23) , C 2 and v are constants . The dependence of 

the fractional species A on time and temperature is clearly 
displayed by this equation . An examination of the data on the MFe 
center in InP validates the above analysis. This can be seen from 
an examination of figures 6.12 and 6.13 which chart out the peak 
heights in the DLTS spectra of the A species at two different 
temperatures (T =175K and 195K). 

6-6: Conclusion 

In this section we first address ourselves to two issues : ( i ) 
What are the experimental signatures in the capacitance transient 
spectroscopy of a broadened gap level as opposed to a sharp delta 
function like gap level ? (ii) To what extent would a standard 
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DLTS analysis of a broadened defect level be in error ? 

Regarding the first question, there are at least three 
indicators of the broadened gap level. First, if the system under 
study is patently inhomogeneous, then the qualitative arguments 
presented in section 6.1 indicate that the level should be 
broadened. Semiconductor alloys are typical suspect systems. An 
extreme example is GaAsSb since it is a semiconductor alloy with a 
large lattice constant mismatch between GaAs and GaSb. Amorphous 
semiconductors such as a-Si:H would be another set of suspect 
systems. Second, an examination of the capacitance transient is 
recommended. If it is non-exponential then one must presume 
broadening. In fact, it was this trait which led earlier workers 
to posit broadening* What we have demonstrated in this work is 
an analytic proof of a long-lived non-exponential character of the 
capacitance transient both for weak disorder (equation (6.7) and 
figure 6.1) and strong disorder (equation (6.12) and (6.13) and 
figure 6.3). Third the DLTS signal is itself broadened over and 
above its normal value. Murawala et al*" have proposed a 
phenomenological relationship connecting AT/Tm the full width at 

half maxima to 


O' AT/T™ = 0.1 + (y/Eo 

This relation may help one to extract a rough value of c. Lastly, 
the Arrhenius plot provides little clue to the existence of 
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broadening. As figures 6.4 - 6.6, 6.9 and 6.10 indicate, the 

Arrhenius plot is linear over four to five time decades, even in 
the presence of significant broadening. 

Regarding the second question, the following conclusions may 
be gleaned from our work. 

(i) The conventional DtTS analysis overestimates the 
activation energy for severely broadened levels. This 
overestimation is often quite small. It may, however, lead to 
small but systematic discrepancies between DLTS and other 
experimental probes such as optical studies. Hall effect etc. In 
the case of a-Si:H studied in section 6.4, the input value was E© 
= 0.82eV while the extracted value was E©' = 0.87eV. The precise 
factor responsible for this discrepancy can be traced to the 
expression for the time constant 

1 _ ln(t2/ti) F(o-/kTni.<?/E©) 

Tm (tZ -tl) 

F<^/kTm,<y/Eo) < 1. 

The functional form of F for small disorder is given by 
(6.9). For large disorder, the complicated form of the 
(6.12) prohibits a simple derivation. F, nevertheless, is always 
less than unity. This character of F is responsible for the 


(6.24) 


(6.25) 


equation 

equation 


over estimation . 
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(ii) The error involved in a standard DLTS analysis based on 
equation (6.3) is minor if (a) o/WYm < i and o-ZEo < 0.1, (b) cy/kTm 
> 1 but cy/Eo < 0.1 and A > 10 s K Thus, standard DLTS analysis 
is fairly robust with respect to disorder. 

We would like to end with a few speculative remarks. The 
power law decay form in equation (6.13) was derived under the 
assumption of additive contributions embodied in equation (6.2). 
This implies that the carriers are emitted independently into the 
relevant band. Such a process is deemed parallel relaxation. 
Mutual correlations and subsequent modifications of the defect 
density of states during the emission process is ignored. 
Preliminary investigation in the relaxation of glassy phenomena 
indicate that non-exponential relaxation arises as a natural 
consequence of serial as opposed to parrallel relaxation 
Similar considerations in the study of capacitance transients may 
yield interesting results. The form of the defect density of 
states has been assumed to be gaussian (equation (6.3)). This form 
is premised on the central limit theorem. A microscopic theory for 
the defect density of states in semiconductor alloys should be 
pursued Finally, detailed analytical and numerical exercise 
similar to the current one must be carried out for other 
experimentally studied quantities of interest (e.g. Hall effect, 
thermally stimulated capacitance etc) to ensure that the 
conventional interpretation is valid for tne disorder broadened 
case. The results of this chapter has been reported in references 


11, 17 and 18. 
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Appendix 

We will indicate the derivation of the modified DLTS 
expression (6.9). Differentiating equation (6.7) 


1 a C(t) 
T a T 


If 


exp [-t/To + 0 . 5(o't/kTTo)^] x 

{[t/To - (o't/kT'ro)^]aTo/aT - (o-t/kTTo)VT} 

I aro/at I >> to/t 


(A.l) 

(A. 2) 


then the last term on the RHS 'of (A.l) may be neglected. This 
condition translates to 


2 + Eo/kT >> 1. 


(A. 3) 


Such a condition is easily met. 

We then exponentiate the term in the curly brackets of (A.l) 
and expand the subsequent logarithm: 

t/To - (c/t/kTro)^ = t exp Cln( l-(o'/kT)^t/To]/To 

% t exp[-(cy/kT)^t/To] (A. 4) 

We obtain next the DLTS extrema by differentiating equation (6.8) 
and using equations (A.l) and (A. 4) for the signal. This 
expression is trivially equation (6.9). 



6-23 


REFERENCES 

1. D.V.Lang, J. Appl. Phys.15.3014 - 3023 (1974). 

2. P.M. Mooney, N.S. Caswell, P.H. Solomon and S.L. Wright, Materials 

Research Society Meeting, San Francisco 15 -18 April (1985). 
T.N.Theis, T.F.Kuech, L . T . Palmateer and P.M. Mooney, Gallium 
Arsenide and Related Compounds, Inst. Phys.Conf. Ser. 74, Ch.4, 

p 241 (1984). 

3. P.Parayanthal and F.H. Poliak, Phys. Rev. Lett. 1822 (1984). 

4. L.Samuelson, S. Nilsson and H .G.Grimmeiss, Phys. Rev. Lett. S3. 
1501 (1984). 

5. P.Omling, L.Samuelson and H . G . Grimmeiss J. Appl. Phys. SI. 5117 
(1984). 

6. P . A . Murawala, Vijay A. Singh, S . Subramanian , S . S .Chandvanker 
and B.M.Arora, Phys. Rev. B 23. 4807 (1984);Proc. 13*’'^ Int. Conf. 
on defects in semiconductors, ed . L.C. Kimerling and J.M.Parsey Jr 
(New York: AIME) p 1165 (1985). 

7. Vijay A. Singh (unpublished). 

8. M.Kanieswka and J.Kanieswka, Solid State Commun. S2. 485 (1985). 

9. D.V.Lang, J.D. Cohen and J . P . Harbison , Phys. Rev. B 23. 109 

(1982). J.D. Cohen, D.V.Lang, J.P. Harbison and A.H.Seregent Solar 
Cells la 119 (1983). 

10. M.Bleicher and E. Lange Solid State Electron. IS. 375 (1973). 

11. Raj K. Singh, Vijay A Singh, J . W .'Corbett , and Amita Das, 

J.Phys.C la ,2177 (1986). 



6-24 


12. C.M. Bender and S.A.Orszag (ed.) Advanced Mathematical Methods 
for Scientists and Engineers (New York; McGraw-Hill) p 267 (1978). 

13. Raj K. Singh, Vijay A. Singh and J.W. Corbett, Semiconductor 
Science and Technology, Z, 716-725 (1987); Z, 726-731 (1987). 

14. M. Levinson, J. Appl. Phys 2638 (1985); M. Levinson 

(private communication). 

15. R.J. Palmer, D.L. Stein, E. Abrahams and P . W . Anderson , Phys. Rev. 
Lett. aZ 958 (1984). 

16. V. A. Singh in Elestronic Structure and its applications. 
Lecture notes in Physics, Springer Verlag, ZSZ, 201 (1987). 

17. Amita Das , Vijay A. Singh , and D.V.Lang, Semicond . Sci . and 

Tech., a 1177-1183 (1988). 

18. "Transient Spectroscopy and Disorder", Vijay A. Singh and Amita 
Das, Radiation Effects and Defects in Solids, (1989). 



6-25 


FIGURE CAPTIONS 

Figure 6.1: Plot of the capacitance transient for small o-/kT (ss 
0.2) and t = 3.4418 s . Note the slower decay of the broadened 
gap level as opposed to the exponentially dictated form. T = 

180K. , exact (equation (6.2)); +++ , aprroximate (equation 

(6.7)); .exponential (equation (6.6)). 

Figure 6.2: Simulated DLTS signal for the parameters of figure 
6. land rate window tz/ti = 4. The differences between the various 
expressions are very small, thus vindicating the standard DLTS 

analysis based on equation (6.3). , exact( equation (6.2)); +++, 

approximate (equation (6.7)); , exponential (equation (6.6)). 

Figure 6.3: Plot of the capacitance transient for = 0.85eV, A = 
4 x10*s”*K~^. The broadening factor o - 0.175eV and T = 340K, 

leading to o'/kT = 6. The transient is long lived and distinctly 
non-exponential. The approximate expression (6.12) is in 
reasonable agreement with the exact one. The parameters selected 

for this plot are similar to ones encountered in a-Si:H . , 

exact (equation(6.2)); +++, approximate (equation (6.12)); 
exponential ( equation( 6 . 6 ) ) . 

Figure 6.4: The Arrhenius plot for a simulated DLTS signal with 
rate window b^/t^ = 4. A typical broadened deep level was 

considered and the parameters employed for the simulation of 
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equation (6.2) were = O.SeV, A = and = O.leV. The 
Arrhenius plot is linear over four time decades (0.1 ms to 1 s). 
Note the difference between the extracted parameters (E ,A). E ' = 

O O 

O.SlSeV, A' = l.SxloW^. 

Figure 6.5: The Arrhenius plot for a typical broadened deep level 
with broadening cf ~ 0,25eV. The plot is linear over four time 
decades (0.1 ms to 1 s). Note the difference between the extracted 
activation energy E ' = 0.607eV, and E = 0.5eV. A' = 2.5 x 

O O 

Q — ^ — >2 

10 s K . Other parameters employed are the same as figure 6.4. 

Figure 6.6: The Arrhenius plot for a typical broadened deep level 
with broadening cr = 0.4eV. The plot is linear over four time 
decades (0.1 ms to Is). Note the difference between the extracted 
activation energy E^' = 0.72eV, and E^ = O.SeV. A' = 2.4 x 
10*^s~*K~^. Other parameters employed are the same as figure 6.4. 

Figure 6.7: Plot of the extracted activation energy E^' against 
the dimensionless broadening o'/E^. The parameters used are the 
same as figures 6.4 - 6.6. There is a systematic overestimation of 
the activation energy for cr/E^ < 0.05. The inset displays E^' in 
the weak disorder limit (cr/E^ > 0.1). 

Figure 6.8: Comparison of the experimental DLTS signal with the 
simulated one for a-Si:H. The parameters are identical to figure 
6.3. The rate window t^/t^ = 4. , experiment; , simulated. 
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Figure 6.9: The Arrhenius plot for a-Si:H. The parameters are the 
same as in figure 6.3 and 6.8. The plot is linear over four time 
decades (0.1 ms to Is). There is an overestimation of the 
extracted activation energy :E ' = 0.87eV while E = 0.82eV. A' = 

o o 

4.4 X 10®s“V^. 

Figure 6.10: The Arrhenius plot for an E3-like level in GaAs Sb 

1 — X X 

(x = 0.14). The plot is linear over four time decades (0.1 ms to 
Is). There is a negligible difference between the extracted 
activation energy E ' = 0.25eV and E = 0.23eV. A' = 3.4 x 10* 

O O 

s~\~^ . 

Figure 6.11: (2a) Generalized configuration coordinate diagram for 
a defect with two metastable states A and B . The shaded region 
around the coordinates Q and Q_ indicates the disorder induced 

A. B 

broadening to a finite DOS from a sharp delta function like defect 
level . (2b) The forward (V^) and reverse (V ) bias during a DLTS 
experiment . The transformations of the defect between the A and B 
configurations are indicated. 

Figure 6.12: The peak height in the DLTS spectra of the A species 

for the MFe center as a function of the logarithm of the filling 

pulse t for T = 175K .The experimental points are taken from 
P 


reference 14. 
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Figure 6.13:The peak height in the DLTS spectra of the 
for the MFe center as a function of the logarithm of 
pulse tp for T = 195K .The experimental points are 


A species 
the filling 
taken from 


reference 14. 
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CHAPTER 7 : CONCLUSION 

In this thesis we have focussed on two issues involving 
disordered semiconductors. These are: (i) the presence of a gap 
level due to a point defect in a disordered medium and (ii) the 
effect of disorder on the analysis of experiments on the basis of 
expressions valid and derived for the ordered case, as is the 
usual practice. 

In chapter 1 we have laid down the broad conceptual framework 
in which the thesis work was carried out. The Coherent Potential 
Approximation^ (CPA) was described using a few well chosen 
examples^. The Koster-Slater^ approach to a gap level was 
illustrated with a toy model for a semiconductor. 

The problem of predicting a gap level in an alloy 
semiconductor may be partitioned as follows: (i) Determining the 
electronic structure of the host which is substitutionally 
disordered. The sophisticated CPA was employed to achieve this 
objective. (ii) Locating the gap level by recourse to a 
Koster-Slater equation. In chapter 2 we have illustrated the 
procedure with the single band model"*. In chapter 3 we first 
develop the CPA formalism for a realistic multiband semiconductor. 
We next detail the application of the Knster-Slater equation. We 
have shown analytically as well as numerically that the prediction 
of the trap level by CPA and the simplest mean field theory VGA 
(Virtual Crystal Approximation) are identical for the vacancy. 
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This suggests that a perturbation approach to strong defect 
potentials may be gainfully employed. We call this the Deep Level 
Approximation (DLA). We have pointed out that since the VGA is 
basically a weak scattering theory, it agrees with the CPA only in 
the domain of weak scattering. 

We have then concentrated on an archetypical defect in 
Ga^_^Al^As: the DX center . This has been widely discussed on 
account of its puzzling physical and chemical properties. We have 
reviewed some of its salient features. It is believed to be a 
donor related level in the alloy Ga, A1 As. The shallow to deep 
instability of this particular defect when the A1 concentration is 
increased beyond 20% led us to view the level as a simple donor 
which gets deepened due to disorder in the environment. We have 
self consistently accounted for disorder using CPA**. Our results 
for a simple substitutional donor in the disordered environment of 
Ga A1 As agree well with the Hall effect measurements. The 

l-X X 

general trend of the level with the A1 concentration 'x' also 
compares well with the experimental value. We have calculated the 
defect level for most of the group IV and group VI donors and have 
compared the theoretical values with those obtained from the 
experiments. Furthermore, our calculations reveal that the defect 
possesses an a^ symmetry, which explains the phenomenon of 
persistent photoconductivity (ppc) associated with this defect. 

The DX center has been identified in pure GaAs under 
pressures exceeding 20 kbar. We conjecture that when pressure is 
applied to a system, it results in a distortion of the lattice. 
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This is then the case of positional disorder. Thus, deepening of 
the level due to pressure can once again be attributed to 
disorder. This conjecture however needs to be explored. A cluster 
type calculation may be gainfully applied to this end. 

The EPR invisibility of the DX center has led some authors to 
conjecture that this defect is a negatively charged donor. We have 
explored this viewpoint using Effective Mass Theory for helium 
like impurities^. The results of our calculations do not rule out 
such a possibility. 

We would like to conclude with the following remarks for the 
DX center. Our calculations^ reveal that disorder can result in 
deepening a level. We have shown it specifically for simple donor 
levels in Ga A1 As. The negatively charged donor, which explains 

1~X X 

EPR invisibility is not ruled out by an Effective Mass type 
approach. Similar effects of deepening of the level due to 
disorder in the host environmfent can be observed on negatively 
charged donors also. However, a firm statement in favour of the 
band model for the DX center cannot be made. This is because the 
cavity radius remains an empirical parameter in our EMT 
calculations. A simple self consistent procedure to delineate the 
dielectric screening £(r) needs to be introduced. Further, the 
presence of slight lattice distortion can also explain the 
phenomenon of divergent optical and thermal activation energies. 

The Effective Mass formalism for a helium like impurity has 
been extended to several other problems'^. We have studied the 
chalcogen impurities in silicon, the anti-site defect in GaAs 
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(which is currently viewed as a popular model for the EL2 center 
in this semiconductor). Our results are in good agreement with the 
experimental values. 

Finally, in chapter 6 we have concentrated on the analysis of 
some experiments. We have demonstrated analytically that a non 
exponential decay in a capacitance transient could be an 
experimental signature for a broadened gap level. We have further 
shown that in the conventional analysis of some experiments like 
DLTS and PATS, the expression employed have been derived for a 
sharp level. This results in an overestimation of the activation 
energy for severely broadened levels®”'*'^ . This may result in a 
systematic discrepancy between DLTS, PATS and other experimental 
probes. It could also result in serious disagreements between 
experiment and theory. 

We would like to conclude that an attempt must be made to go 
beyond the Koster-Slater paradigm. Since several centers in the 
alloy semiconductors are metastable, total energy calculations are 
warranted. We have taken a step in this direction in chapter 2, 
where we use the Hellmann-Feynman theorem to calculate the total 
electronic energy for a single band tight binding Hamiltonian*. 

On the other hand an attempt ought to he made to go beyond 
the mean field descriptions of the disordered host. The Coherent 
Potential Approximation is a powerful approach. But localization 
effects are largely washed out by it. It would be a moot point to 
explore the interaction of a 'localized' defect wave function with 
the host eigenfunctions which are not extended Bloch states but 
are Anderson^^ localized. 
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